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Abstract 



In this PhD thesis we analyze the medium modification of the iroop coupling and 
subsequently discuss its possible influence on the dilepton production rate in relativis- 
tic heavy-ion collisions. 

The first part of the thesis is devoted to the medium modifications of the uo — > 
7r°7*, p a — > 7r a 7* and 7r° — > 77* decays in nuclear matter. We use the relativistic 
field theory formalism describing modifications of the coupling by one-loop Feynman 
diagrams with nucleons and A(1232) isobars. For simplicity we work in the leading 
baryon density approximation and at zero temperature. We consider cases when the 
decaying particle is at rest or at motion with respect to the medium. The kinematics 
is considered in the rest frame of the medium. In the medium, the structure of the 
Troop amplitude is more complex than in the vacuum, which is due to the presence of 
the additional four-vector describing the velocity of the medium. We investigate the 
dependence of medium effects on the dilepton invariant mass and observe a sizeable 
increase of the effective coupling constants compared to their vacuum values. The 
-nujp coupling constant increases by about a factor of 2 for the u — > 7r°7* decay and by 
about a factor of 5 for the p a — > 7r a 7* decay. The effect grows with density. Similar 
conclusions are obtained for a particle moving with respect to the medium. 

In the second part of the thesis we apply our model to evaluate the influence of the 
medium effects on the dilepton production from the Dalitz decays. In the calculations 
of the dilepton spectrum we use the model of the hydrodynamic expansion of the fire 
cylinder, including the longitudinal and transverse expansion. In our comparison to 
the CERES dilepton data we take into account the experimental cuts, which is very 
important in the detailed numerical analysis. The medium modifications enhance the 
effect of the Dalitz decays, however, the numbers are still below the experimental 
data. 

It is worth to point out that in the medium the p decay due to the isospin degen- 
eracy becomes as important as the uj decay. The main new and original results of the 
thesis are: 

1. Significant in-medium increase of the nup coupling constant. 

2. Large enhancement of the Dalitz uj — > 7r°7* and p a — > 7r a 7* decays in the range 
0.2 — 0.6 GeV of the dilepton invariant mass. 



Acknowledgements 



First and foremost I have the pleasure to thank my supervisor Dr. Wojciech 
Broniowski for guiding me through my PhD. 

I am also very grateful to Prof. Wojciech Florkowski for many valuable comments. 

I would like to thank Prof. Edward Kapuscik, the head of International PhD 
studies at IFJ PAN for help and for many interesting discussions. 

Many warm thanks to the staff and colleagues of Theory Department in the In- 
stitute of Nuclear Physics for numerous discussions during my PhD studies. 

I warmly thank my sister Marta and Steven Steinke from University of Arizona 
for linguistic help. 

Finally I would like to thank my family for support, and all people who in some 
way have helped and encouraged me. 
Special thanks to Piotr Czerski. 



Contents 



Introduction 



I Modification of the irup coupling in nuclear medium 



1 


Modification of Hadrons in medium 






1.1 The Walecka model 








1.2 Ouark condensate 


in medium 






1.3 Brown- Rho scaling 






1.4 OCD sum ruled 




1.5 Predictions from disoersion relatione 






1.6 Lowering; of the Dosition of o° peak! 




1.7 Modification of coupling constants! 



2 


TTuio vertex in nuclear medium 








2.1 Our model 








2.2 Decavs of Darticles at rest with resDect to the medium 






2.3 Decav of particles moving with respect to the medium 



3 Summary of Part I 



II Dilepton production rate 



4 Electromagnetic sienals from hot and dense matter 


4.1 Overview of experimental measurements 




4.1.1 CERES /NA45 experiments 


4.1.2 HELIOS experiment] 


4.1.3 NA38/NA50 experiment 


4.1.4 Dilepton spectra at BEVALAC energies 


4.1.5 Other experiments 



9 



4.2 Overview of theoretical model! 



5 


Formalism of the dileDton Dro 


duction from vector mesons 






5.1 Vector Dominance Model (VDM) 










5.2 Dalitz decavsl 




5.3 Model of the hvdrodvnamic exoansion of the fire-cvlinder 






5.4 Kinematics for dilepton vield 





65 

69 

70 



6 Results of our model 



7 Summary and discussion 



8 


ADDendices 






8.1 ADDendix A 






8.2 ADDendix B 






8.3 ADDendix C 






8.4 Appendix D 





83 
87 

91 

91 
92 
94 
96 



10 



Introduction 



Physics of ultra-relativistic heavy-ion collisions is a quickly developing branch of sci- 
ence which allow us to investigate the properties of nuclear matter in terms of ele- 
mentary interactions. The production of dileptons (e + e~, p + p~) is one of the most 
important probes in studying the dynamical evolution of the nuclear collision pro- 
cesses. Since they do not interact strongly, the dileptons escape unthermalized from 
the hot and dense matter during all stages of the evolution. A major source of dilep- 
tons comes from the direct and the Dalitz decays of meson resonances in particular 
the u, p, 7], rf or 0. The properties of vector mesons such as masses or widths, are 
known to be significantly modified in nuclear matter, PJ El EJ- These in- medium 
properties play an important role in the low mass dilepton production region, namely 
between 0.2 — 0.6 GeV of invariant masses, covered recent heavy-ion collision exper- 
iments, |H El IE]- The dilepton production enhancement in central nucleus-nucleus 
collisions in the low-mass region has been a topic of great interest in the last years, 
[3 IH1 M EH IU| • This phenomenon is frequently explained by assuming that the vec- 
tor meson spectral functions undergo substantial modifications through the strong 
interactions with the nucleons in the medium. Many research groups search for a 
theoretical description of the low-mass dilepton enhancement however, a convincing 
explanation has not been reached and the problem still remains a major puzzle in 
relativistic heavy-ion physics. 

Trying to explain this relevant problem, in this thesis we investigate modifications 
of the nup coupling constant in the context of the uj and p Dalitz decays. We analyze 
these decays in nuclear matter applying a relativistic hadronic framework incorpo- 
rating nucleons and A (1232) isobars. We find that the nujp coupling constant for 
uj — > 7T°e + e~ and p a — > ir a e + e~ decays are considerably enhanced by the medium as 
compared to the vacuum values. We observe an enhancement by about a factor of 2 
for the uj meson and by about a factor of 5 for the p meson. We have to stress that the 
Dalitz decay of the p meson is found to be equally important to the omega decay, ^2j . 
The in-medium-modified ixujp coupling constant increases the corresponding widths 
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of considered decays. In order to estimate the dilepton yields we use the vector domi- 
nance model (VDM) and include the effect of the expansion of the medium. Next, we 
compare our model results to the dilepton yield from the CERES Pb+Au experiment 
in the low-mass region carefully including the experimental kinematic cuts. Finally, 
our results show that medium effects from the u and p Dalitz decays are significant 
and about two times larger compared to the vacuum. 

Highlights of the experimental situation 

The physics of ultra-relativistic heavy ion collisions, i.e., collisions of atomic nuclei 
in which the center-of-mass energy per nucleon is much larger than the nucleon rest 
mass, has been studied experimentally with immense intensity. First experiments have 
been performed at the Alternating Gradient Synchrotron (AGS) in the Brookhaven 
National Laboratory (BNL) and at the Super Proton Synchrotron (SPS) at CERN, 
where ion beams were accelerated with the center of mass energies around y/s ~ 
5 A GeV and y/s ~ 20 A GeV, respectively. 

In 2001 the first data were collected from the Relativistic Heavy Ion Collider 
(RHIC) at BNL at the energy of y/s ~ 200 A GeV. Between the years 2001-2004 the 
next three runs took place. In the near future heavy ions will also be injected into the 
CERN Large Hadron Collider (LHC) reaching the energy of y/s = 10 A TeV. The 
principal goal of experimental and theoretical initiative in ultra-relativistic heavy-ion 
collisions is the observation and understanding the physics of hot and dense medium. 

Many probes have been proposed to analyze the behavior of hot and dense hadronic 
matter. Among these probes the dileptons stand out, because they couple directly 
to vector mesons and interact only electromagnetically. Therefore they are likely to 
bring information about the hot and dense matter formed in the early stages of nu- 
clear collisions. A schematic view of characteristic dilepton sources in ultrarelativistic 
heavy-ion collisions is given in Fig. ^ One can distinguish three basic regions with 
respect to the invariant mass of the dilepton pair. 

In the high-mass region, M e+e _ = 3 GeV, dileptons come dominantly from hard 
processes such as the Drell-Yan annihilation occurring in the early stage, where the 
excited hadronic system is far from thermal equilibrium. After rapid thermalization 
the QGP phase may be established, where the dilepton production proceeds predom- 
inantly via the quark-antiquark annihilation. Upon expansion and cooling, the QGP 
is converted into a hot hadron gas. At this stage dileptons are radiated from pion 
and kaon annihilation processes as well as other collisions between various hadrons. 
In contrast to the light vector mesons (u, p), the lifetime for the heavy quarkonium 
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Figure 1: Major dilepton sources as a function of invariant mass in ultrarelativistic 
heavy-ion collisions. In the high-mass region above invariant masses M e+e _ = 3 GeV, 
dileptons come mostly from hard Drell-Yan annihilation processes. Here strong J/\& 
suppression and abundance are interpreted as a sign of color deconfinement. In the 
intermediate-mass region between 1 — 2.5 GeV, the thermal signal from QGP could 
be related to the observation of associated DD production. In the low-mass region 
below 1 GeV, dileptons originate from Dalitz decays of neutral mesons such as it , rj, 
r]' — > e + e~7, uj — > 7r°e + e~, and from direct decays, u, p, <fi — > e + e~. This region is 
very sensitive to hadron modifications in the nuclear matter. 

states such as J/\I/ and T is substantially longer than the typical life-time in the 
hadronic fireball. Therefore these mesons will mostly decay after freezeout, and hence 
will not feel the medium. 

The intermediate mass region, above ~ 1.5 GeV, might be proper to observe a 
thermal signal from the plasma radiation, because QGP can only be formed at higher 
temperatures than a hadronic gas. This signal could be revealed by the observation 
of the associated DD production in this intermediate region. 

Next, when the freezeout stage is reached, the dominant sources are hadronic 
resonances as well as Dalitz decays, mostly from ir°, rj, u, and p mesons. The two- 
body annihilation processes are dynamically enhanced through the formation of light 
vector meson resonances such as the p, u, and mesons, which directly couple to 
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lepton pairs. Thus, the invariant mass distribution of the lepton pairs reflects the mass 
distribution of the vector mesons at the moment of decay. Therefore, vector mesons, 
possibly modified by medium effects, play an important role in the measurement of 
dilepton rates in relativistic heavy-ion reactions. All these above-mentioned dominant 
sources feed into the low mass region, Mz + /_ < 1 GeV, which is particularly sensitive 
to in-medium modifications of the light hadrons. 

The experimental measurements of the dilepton spectra in relativistic heavy-ion 
collisions have mainly been carried out at the CERN SPS by the CERES/NA45 
(low-mass region) |4*1 [T3"l IT4*j . HELIOS-3 (low- and intermediate-mass region) 5J and 
NA38 /NA50 (intermediate- and high-mass region) ^SJ CZ] • In the near future, at 
RHIC and CERN-SPS, dilepton spectra will be measured by the PHENIX and NA60 
[T5] collaborations, respectively. At much lower bombarding energies, dilepton data 
have also been taken by the DLS collaboration at SIS (GSI) f01l2DI 5 where only the 
low-mass region up to the kinematical limit of around 1 GeV may accessible. These 
experiments focused on the role of high baryon density. 

Recently, the CERES jl] and HELIOS 5J collaborations have observed that cen- 
tral nucleus-nucleus collisions exhibit a strong enhancement of low-mass dilepton pro- 
duction as compared to proton-nucleus reactions. The data analysis is based on the 
so-called hadronic cocktail model which works very well to describe the proton-nucleus 
data. On the other hand, the hadronic cocktail strongly underestimates the dilepton 
yield from nucleus-nucleus collisions. The model is based on data from proton-induced 
collisions which are scaled with the multiplicity of charged particles. A lot of efforts 
have been undertaken in order to understand phenomenon of the low-mass enhance- 
ment. 

In theoretical attempts to address the problem, the increased dilepton yield in 
nucleus-nucleus collisions was attributed to enhanced p-meson production via 7r + 7r~ 
annihilation. Thus, many theoretical groups, j2H 122 I2H| have included the p — > 
7r + 7r~ process within different models for the space-time evolution of nucleus-nucleus 
reactions. Their calculations are based on the vacuum meson properties, and all 
results are in disagreement with the data. This has led to the suggestion of various 
medium effects, such as the reduction of the vector meson masses, that might be 
responsible for the observed enhancement. Brown and Rho suggested the hypothesis 
that the vector meson masses drop in the medium according to a simple scaling law, 

mm 

There are also many approaches based on QCD sum rules OEH], which confirm 
the scaling hypothesis. 
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The conclusion from the existing studies is that although conventional mecha- 
nisms are sufficient to account for the dilepton spectra from proton-induced reactions, 
medium effects are needed to explain the low-mass dilepton enhancement observed 
in nucleus-nucleus reactions. Thus, in this thesis we examine the in-medium nup 
vertex and its large influence on the dilepton production rate in relativistic heavy-ion 
collisions. We believe that our model together with other medium effects may resolve 
the dilepton puzzle. This is the main reason of our study. 

Organization of the Thesis 

The first part of this dissertation is devoted to the modification of hadron properties, 
in particular the irup vertex, in the presence of nuclear medium. 

As the introduction to the subject, in Chapter 2 we first review some theoreti- 
cal approaches describing the medium modifications of hadron properties. Next, we 
present our original analysis, describing our model (section 2.1). The model is based 
on fully relativistic hadronic theory where mesons and nucleons interact with each 
other and with the A isobars. The presentation concentrates on Dalitz decays such 
as 7T° — > 77*, uo — > 7r°7* and p a — > 7r a 7*. In such processes the inup coupling constant 
appears. We present propagator, vertices for the interactions with and without A 
resonance, and the in-medium structure of the nup amplitude. This structure is more 
complicated because of the additional four-vector, the four-velocity of the medium. 
We analyze the Dalitz decays at rest, section 2.2, as well at motion, section 2.3, with 
respect to the medium. In the case when the decaying particles are moving with 
respect to the medium we consider the kinematics in a convenient and non-textbook 
way. The kinematics is usually analyzed in the rest frame of the decaying particle, 
whereas we consider it in the rest frame of the nuclear medium (with the four- velocity 
of the medium u = (1,0,0,0)). We investigate the dependence of medium effects in 
the dilepton invariant mass and observe their influence on the effective nup coupling 
constant compared to its vacuum value. We point out that the nup coupling constant 
is considerably enhanced in the presence of the nuclear medium for both the uj and 
p Dalitz decays. The uj and p can have two polarizations. When the decaying parti- 
cle moves with respect to the medium each polarization behaves differently, thus the 
properties of decaying particle are different. We calculate the dependence of the width 
for both the transverse and longitudinal polarizations as a function of the invariant 
mass of the dileptons, and show significant medium effects. Finally, in Chapter 3 we 
discuss how nuclear matter influences the ixujp coupling constant in the previously 
mentioned Dalitz processes. 
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In the second part of the thesis, we focus on applying our model to the production 
of dileptons in relativistic heavy-ion collisions. We begin this part with a review of 
major dilepton experiments and some theoretical models connected to the low-mass 
region. This part starts out by describing the experimental measurements of the dilep- 
ton yield. Here we concentrate on the experimental background and discuss relevant 
dilepton experiments, mainly at CERN SPS in the proton-nucleus and nucleus-nucleus 
reactions. We discuss results at high bombarding energies of 158 A GeV or 200 A GeV 
and at the lower energies of 40 A GeV and 1 — 5 A GeV. We focus on the low-mass 
region, because here in particular the experimental results are in disagreement with 
many theoretical models. The experimental data are usually compared to the hadronic 
cocktail model, therefore we present a short description of this approach. In section 
4.2 of Chapter 4, we describe the most popular theoretical explanations related to 
the low-mass dilepton enhancement, based on the in- medium modifications of vector 
mesons. 

In order to estimate the dilepton production rate from the Dalitz decays of vector 
mesons we use the vector dominance model (VDM) reviewed in section 5.1 of Chapter 
5. Section 5.2 discusses in detail the structure of the three-body Dalitz decays. In 
section 5.3 we describe briefly the model considered for the hydrodynamic evolution 
of the fire cylinder from Ref. (2^1 E2j. Next, in section 5.4, we present details of our 
formalism used to calculate the dilepton production from the Dalitz and direct decays 
of the u and p mesons. We use the experimental acceptance cuts in our analysis. It is 
quite nontrivial how to construct and include these cuts in the calculations, therefore 
in many other theoretical works the acceptance cuts are frequently not taken into 
account. 

In section 5.5 we present our numerical results, starting from calculations without 
expansion of the fireball and without the CERES acceptance function. Our plots are 
done for two typical temperatures of the fireball. Next, we present the dilepton yield 
from direct and Dalitz decays with expansion and acceptance cuts. We compare our 
numerical results to the CERES experimental data also for two different temperatures. 

In the medium we observe a significant enhancement for the Dalitz decays of the 
<jj and p mesons in the dilepton mass region 0.2 — 0.6 GeV compared to the vacuum 
case. 

Although the modification of the irup vertex increases the contribution from the 
Dalitz decays, the obtained yields alone fall short of describing the data. Yet, they 
help by contributing more to the theoretical curves, and together with other effects 
not considered here (such as dropping masses, broadening widths, or Dalitz decays 
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of i] and rf mesons), may help to explain the long standing problem of the low-mass 
dilepton enhancement seen in relativistic heavy-ion collisions. 

The new and original material of this thesis is contained in Chapter 2, section 5.4 
and Chapter 6, while Chapter 1, Chapter 4 and sections 5.1, 5.2, 5.3 contain useful 
reviews necessary for the integrity of the text. 
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Part I 

Modification of the iruup coupling in 

nuclear medium 



19 



Chapter 1 

Modification of Hadrons in medium 



As mentioned in the Introduction, in heavy-ion collisions the particle production is 
related to the evolution of hot and dense matter. Heavy-ion collisions provide the only 
way to compress and heat up nuclear matter in the laboratory conditions. From these 
heavy-ion data we obtain important information useful in construction of models for 
the early universe, supernova explosions or neutron stars. In the so-called relativistic 
regime, at energies above 1 GeV per nucleon, heavy-ion collisions are likely to yield 
particle densities of the order of a few times the normal nuclear matter density. In such 
system the particles cannot propagate completely freely (the Compton wavelength is 
comparable with mean free paths of particles). In this situation we expect that in hot 
and dense environment the particle properties, i.e. their masses, widths, or coupling 
constants, are changed and these in-medium modifications may lead to experimentally 
observed phenomena. 

The problem of how the properties of hadrons change in hadronic or nuclear matter 
in comparison to their free values has attracted a lot of attention [21 El El El EH EH 
EaE3EaEDE3E3ElE3E3E3EaE3EailHll2!. Till a few years ago there were 
only theoretical predictions [^HEHIED!, for the in- medium modifications. At present 
we also have some indirect experimental evidence from the dilepton production in 
hadronic collisions, (which we discuss in more detail in the Part II, Chapter 4) j3J E] 
or from data on the 7r + 7r~ correlations from STAR, where the p meson peak moves 
to lower masses [33] • The existing calculations are based on simple scaling of masses 
P [7], numerous hadronic models [21 El El El GUJ E21 GS3 GUI ESI GIHl EH EH EH1 EOl EIJ, 
or QCD sum rule techniques [231E21I32I- More fundamental predictions are based on 
low-density expansion and dispersion relations O ESI ESj • In all these calculations 
properties of hadrons are changed by the presence of the medium through strong 
interactions with nucleons from the medium. 
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Below we will briefly review major theoretical approaches for in-medium properties 
of vector mesons, namely, p, u, and mesons. 



1.1 The Walecka model 

A simple and very illustrative approach showing the effects of medium modifications 
is Quantum Hydrodynamics (QHD) [21 E] in its simplest form known as the Walecka 
model. It is based on the mean-field theory of hadrons, where the scalar and vector 
mean fields generated by the nucleon sources are themselves self-consistently responsi- 
ble for the modification of the nucleon mass. The Lagrangian density for the Walecka 
model is 



where \1/ is the baryon field of mass M, is the neutral scalar-meson field with mass 
m s , is the neutral vector meson field of mass m v , and F^ u = d^V u — d v V^. 
The term SC contains renormalization counterterms. The nucleons interact via the 
exchange of isoscalar mesons with the coupling of the scalar field to the baryon scalar 
density and the vector field V M to the conserved baryon current ^7^^/, obtained 
through minimal substitution. Because the exact solutions to the field equations are 
very complicated, a mean-field approximation is used. In a mean-field approximation, 
the meson field operators are replaced by their expectations values, which are classical 
fields: 



Now the mean-field equations can be solved exactly with the solution becoming in- 
creasingly valid with increasing baryon density. Usually, the mean-field equations are 
solved in two approximations. In the mean-field theory one calculates a baryon self- 
energy which is generated by the presence of all the nucleons in the occupied Fermi 
sea, and the effect of the Dirac sea is neglected. In contrast, in the relativistic Hartree 
approximation one includes the contribution to the baryon self-energy arising from 
the occupied Fermi sea as well as from the full Dirac sea. In consequence, the baryon 
self-energy diverges in the relativistic Hartree approximation and must be renormal- 
ized. In both approximations, the mean scalar meson field 0o is responsible for a shift 
of the nucleon mass M* in the nuclear medium relative to its vacuum value M. In 



£= * [ 7M (^-^)-(M-^0)]^ + i(^0^0-m^ 2 ) 
- iF^F^+^mlV^V + SC, 



(1.1) 



-> ((f)) = O , 

V, - (V M > EE S^Vq. 




contrast to the ground-state expectation value of the vector field, which is determined 
by the conserved baryon density, p#: 

V° = %%°*> = ^p B . (1.3) 

The expectation value of the scalar field, and consequently the effective mass of the 
nucleon, is a dynamical quantity that must be determined self-consistently from the 
equations of motion 

2 2 

M* = M- g s (f> — M — -%<M) — M — ^p s , (1.4) 

trig mj, 

where p s = (fy^) is the scalar density of the nucleons. Since p s > 0, from the above 
equation we see that indeed nucleon mass, M*, is lower in medium than in vacuum. 

The model has already been used extensively in calculations of nuclear matter 
and finite nuclei. The saturation of nuclear matter and the strong spin-orbit splitting 
observed in finite nuclei were among the first successes of the model |4~T) 14 7| l4*Hl I49j. 
The model has also been used with considerable success to analyze such diverse topics 
as collective modes in nuclear matter, isoscalar magnetic moments, and electroweak 
and hadronic responses from finite nuclei. 



1.2 Quark condensate in medium 

The authors of Ref. [50J have studied quark and gluon condensates in nuclear matter. 
These condensates are expectation values of local composite operators such as qq and 
G a ^ u G atlv ', where q is an up or down quark field and G a is the gluon field-strength 
tensor. The authors of Ref. [5U] focused on describing the properties of hadrons in 
nuclear matter in terms of in-medium quark and gluon condensates, which are shifted 
from their vacuum values. The condensates are analyzed in the ground state of nuclear 
matter. Simple expressions for these quantities, which are model independent to first 
order in the nucleon density, have been developed. 

As a result, the ratio of the in-medium quark condensate to its vacuum value is 
given by 

(qq) pk = (qv)vac(i - ^yZ), (i-5) 

where (qq) PN and (qq) vac are the medium and vacuum values of the quark conden- 
sate, m-x is the pion mass and f n is the pion decay constant. The quark condensate 
at low densities is related to the nucleon a term, cr^, which can be expressed as 



(Tat = (N\mqq\N), where m is the current mass of light quarks. Eq. 11.51 shows the in- 
medium quark condensate to first order in the nucleon density. There are higher-order 
corrections, see Ref. [HO]; therefore, the model-independent results are valid at suffi- 
ciently low nucleon densities. To first order in the nucleon density, the relationship of 
the in-medium quark condensate to its vacuum value depends on the values of the pion 
mass (138 MeV), the pion decay constant (93 MeV) and the a term (about 45 MeV). 
The value of nucleon density at saturation is = 0.17 fm~ 3 . Depending on the 
precise value of the a term, the quark condensate is reduced considerably at nuclear 
matter saturation density - it is roughly 25 — 50% smaller than the vacuum value, 
thus the shift in the quark condensate is large. This description is complementary 
to the Walecka model from QCD (M ~ (<?<?))• On the other hand, the finite-density 
effects modify the gluon condensate to a much smaller extent than they do for the 
quark condensate. The decrease in the gluon condensate at nuclear matter saturation 
density is around 5%. Hence, it does not affect the particle properties as much. 



1.3 Brown- Rho scaling 

Brown and Rho (BR) have suggested the hypothesis that the hadron masses drop in 
the medium according to a simple scaling law Ij, given by 

= K = < = ft (16) 
m N m p & 

where /„- is a pion decay constant and the medium-modified parameters are indicated 

by asterisks. This hypothesis is based on phenomenological implementation of the 

restoration of chiral symmetry in the framework of an effective theory. 

At least qualitatively, this scaling conjecture is supported by studies based on the 

QCD sum rule approach, as well as by other hadronic models E3 EH M M E21 HHH 

EH EH ESI GUI EE GIHl EH 



1.4 QCD sum rules 

A natural and direct use for in-medium condensates is made in QCD sum rule cal- 
culations of hadronic properties in nuclear matter. This approach can be used to 
predict in-medium spectral properties (effective masses or self-energies) of baryons 
and mesons. 

The QCD sum rule approach aims to understand the physical current-current 
correlation functions in terms of QCD. It is done by relating the observed hadron 



spectrum to the non-perturbative vacuum structure of QCD. Through the use of 
this technique Hatsuda and Lee, [21] have extracted the medium dependence of the 
non-strange vector-meson masses. They have obtained the following result 

1-0.18(^0, (1.7) 
Po 

l_ .15a(^) (1.8) 
Po 

where a is the nucleon strangeness content, a = (N\ss\N)/(N\qq\N), taken to be 
a = 0.17. The p and uo meson masses decrease significantly with density due to the 
strong density dependence of the light quark condensate. Since the strange quark 
condensate does not change much in nuclear media as a result of the small nucleon 
strangeness content, the meson mass shows a weaker density dependence. The fact 
that these masses decrease as density increases has been taken as a confirmation of 
the Brown- Rho scaling pQ. 

The QCD sum rules for vector mesons were reanalyzed in Ref. [SJ. The uncer- 
tainties in various in-medium quark condensates were assessed using the Monte-Carlo 
method analysis. It was found that at normal nuclear matter density m*/m p = 
0.78 ± 0.08, in nice agreement with the results of Refs. [21]. In Ref. (22] a decrease of 
the p meson mass with increasing density like that of Ref. [21] has also been obtained. 
The predictive power of the QCD sum rule approach was addressed in Ref. [23] tak- 
ing into account the possible modifications of the p meson mass and width. In this 
approach the p meson spectral function was parameterized by a simple Breit-Wigner 
form containing two fitting parameters, the pole mass and the width, and the values 
for mass and width are determined by requiring that the QCD sum rules are satisfied. 
It was found that in the case of vanishing width, the p meson mass indeed decreases 
in nuclear matter. On the other hand if the width increases it is then possible to 
satisfy the QCD sum rules with a constant or even an increasing p meson mass. The 
predictive power of the QCD sum rule approach thus depends on how well one can 
constrain the phenomenological p meson spectral function in the nuclear medium. 



and 
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1.5 Predictions from dispersion relations 

The in-medium mass shift and broadening of the width are due to interactions of the 
hadron with the medium. Thus, one can use phenomenological information on this 
interaction to calculate the mass shifts. The authors of |HS] have used the result that 



the mass shift of a particle in medium is related to the forward scattering amplitude, 
f sc , of this particle on the constituents of the medium. 
This gives the mass shift 

Am = -2ir—Ref sc (1.9) 
m 

where m is the vacuum mass of the particle and p denotes the density of the medium. 
The width broadening is given by 

Ar = —k(T, (1.10) 
m 

where k is the particle momentum and ko = 4irlmf sc . Equations 11.91 and 11.101 are 
correct when the particle's momentum k is larger than few hundred MeV, and Ref sc 
fullfils the requirement |i?e/ sc | < d, (where d is the distance between medium con- 
stituents). Additionally the main part of the scattering should proceed through small 
angles, less than 1, and the medium constituents should have some momentum dis- 
tributions, such as Fermi-Dirac or Bose-Einstein distributions, for finite temperatures 
and chemical potentials. In such cases, in the right-hand sides of Eqs. 11.91 and, 11.101 
an averaging over the momentum distribution of constituents must be performed. 
Thanks to this approach one may formulates explicitly the applicability conditions 
discussed above. When gas in thermal equilibrium plays the role of the medium, the 
equivalents of Eqs. 11.91 and. 11.101 can be derived in the framework of the thermal field 
theory. 

Using this approach, authors of Ref. [SH] have estimated p meson mass shift 
and width broadening in the case of p-mesons produced in heavy-ion collisions, with 
following conclusions: For p-mesons produced in high energy heavy ion collisions the 
mass shift is small, but the width broadening is so large that one can hardly observe 
a p-peak in e + e~ or p + p~ mass distributions. In low energy heavy-ion collisions a 
p-peak may be observed in e + e" or p + p~ mass distributions as a broad enhancement 
approximately at the position of p-mass. 

1.6 Lowering of the position of p° peak 

Experimental indications for in-medium physics also come from RHIC [13] , where the 
data obtained by the STAR collaboration have been interpreted as the lowering of 
the position of the p° meson peak, see Fig. 11.11 This is the first direct measurement 
of p(770)° — ► 7T + 7r _ in heavy-ion collisions. The important observation is that the 
p° mass is lower in peripheral Au+Au collisions than in p+p collisions. Apart from 
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Figure 1.1: The tt + ti invariant mass distributions for peripheral Au+Au interactions. 
The p peak is around 700 MeV, much lower than the vacuum value. 



the dynamical interactions with the surrounding matter, interference between various 
7r + 7r~ scattering channels, phase space distortions due to the rescattering of pions 
forming p°, and Bose-Einstein correlations between p° decay daughters and pions 
in the surrounding matter are among possible explanations for the apparent change 
of the p° meson. Hopefully, more detailed information will be provided by future 
correlation experiments. 



1.7 Modification of coupling constants 

Meson properties are changed due to strong interactions with nucleons of the medium. 
If the mass and widths of hadrons can be significantly modified, one can expect 
that also the coupling constants are altered. There are numerous studies of mesonic 
two-point functions in the literature, Ref. [S3], which in general contained medium 
modified mass M* and medium modified width T*, but only a few devoted to meson 
three-point functions. In their studies of the p meson in-medium spectral function in 
Ref. jH], Herrmann, Friman and Norenberg have analyzed the pint vertex for the p 
at rest with respect to the medium. They have considered the above vertex on the 
basis of a hadronic model with the A isobar and with nonrelativistic couplings. Song 
and Koch in Ref. [HE] have taken into account the temperature effects on the pirn 
interaction. Krippa in Ref. [EE] has computed the effects of density on chiral mixing 
of meson three-point functions. The authors of Ref. (EH EHl EDJ ED] have analyzed the 



to — > 7T7T and p — > 7T7T decays in nuclear medium. They have studied the pirn vertex in 
the fully relativistic framework with relativistic interactions, incorporating nucleons 
and A (1232) isobars, and found that the medium effects on the pirn constant are large, 
and dominantly come from the processes where the A is excited in the intermediate 
state. This has immediate consequences for the direct decay of the p meson into 
dileptons (such as broadening of width). 

An important factor brought in by the presence of the medium is that processes 
which are forbidden in the vacuum by symmetry principles are now made possible. 
In this work we concentrate on the nup coupling constant and analyze its significance 
for the Dalitz decays of vector mesons. 



Chapter 2 

nup vertex in nuclear medium 



So far we have discussed major theoretical approaches related to the medium modi- 
fications of particle properties. In this chapter we begin the original material of the 
thesis. We focus on the irup coupling constant and its medium modifications. The 
nup vertex appears in such processes as: u> — > 7T°7*, p a — > 7r a 7*, and n° — > 77*, see 
Fig 12.11 We expect that medium effects may have a significant influence and may 




Figure 2.1: The physical processes where the iruop vertex appears. Wavy lines indicate 
the p orw, dashed lines the pions, the hatched blob indicates the medium-modified 
vertex, the black dot is equal to or — - as follows from the Vector Dominance 
Model (VDM). 

affect the existing calculations based solely on the vacuum values. Our calculation of 
the in-medium irup vertex is made in the framework of a fully relativistic hadronic 
theory, where mesons interact with nucleons and A isobars. Among the baryon reso- 
nances, the A is most important due to the small A — N mass splitting and a large 
value of the 7riVA coupling constant. In our investigation we work at zero tempera- 
ture to obtain the vertices of Fig. 12.11 at least as a first approximation. We use the 
leading-density approximation, which makes the calculation simpler, as no integration 
over the nucleon momenta is necessary. This Chapter is divided into three sections. 
In section 2.1 we present our model with all the needed diagrams, propagators and 
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vertices for interactions of mesons with nucleons and with A resonances which we use 
in our calculations. Also, the tensor structure of the nup vertex is presented. Next, 
we look at how the effective coupling constant is changed compared to the vacuum 
value. First, in section 2.2 we analyze our decays in the rest frame with respect to 
the medium. Later, in section 2.3, keeping non-zero three momenta of the us and p 
we look separately at the longitudinal and transverse polarizations. In both cases we 
observe that the effective coupling constant is significantly enhanced in the presence 
of nuclear matter. 

2.1 Our model 

In our analysis we apply a conventional hadronic model with meson, nucleon and 
A degrees of freedom. This model is similar to the Walecka Model, which we have 
discussed in Chapter I. In Fig. 12.21 the considered diagram for the irujp coupling 
consists of two parts, the vacuum one and the induced by the medium. 




Figure 2.2: The full Troop vertex consisting of two parts, the vacuum part and the 
medium part, respectively. 

There are three types of in-medium diagrams, see Fig. 12.31 diagrams with three 
nucleons, Fig. 12.31 (a), two nucleons and one A isobar, Fig. 12.31 (b), and with one 
nucleon and two A isobars, Fig. 12.31 (c). 

The diagrams of interest are those which include at least one nucleon, because 
only these depend on the baryon density, carried by the Fermi sea. Since the A isobar 
has a large coupling constant it has a large contribution to the considered processes. 
The diagram with three A is not considered, because such a diagram corresponds to 
vacuum polarization effects, which are not considered. 




(a) (b) (c) 

Figure 2.3: In-medium diagrams included in the present calculation (crossed diagrams 
not shown). Solid lines indicate the in-medium nucleon propagator, double lines the 
A propagator. 

The in-medium nucleon propagator 

Next, we introduce the in-medium nucleon propagator and the Rarita-Schwinger A 
propagator which we use in our calculations. We show also all the needed vertices, 
and the in- medium tensor structure of the in- medium nup coupling constant. 

Solid lines in Fig. 12.31 denote the in-medium nucleon propagator, which can be 
conveniently decomposed into the free and density parts [2|: 

iS(k) = iS F (k)+iS D (k)=i(Yki l + m N )[— - 2 , + 

k z — m z N + is 

—6(k -E k )e(k F -\k\)\, (2.1) 

where denotes the nucleon mass, = \Jm 2 N + k 2 , and kp is the Fermi mo- 
mentum. The density effects are produced when exactly one of the nucleon lines in 
each of the diagrams of Fig. 12.31 involves the nucleon density propagator, Sp,. The 
diagrams with more than one Sp> vanish for kinematic reasons, see Appendix D. On 
the other hand, diagrams with only Sp propagators do not contain Fermi-sea effects. 
Therefore, these vacuum diagrams are not considered in the present study. Finally, 
we have only one possibility of having exactly one nucleon density propagator in each 
of the diagrams of Fig. 12.31 

The Rarita-Schwinger propagator for the A baryon 

In relativistic field theory the covariant description of spin 3/2 particles is based on 
the Rarita-Schwinger formalism [HT], where the fundamental object is the spin- vector 
field The most general Lagrangian for the massive spin 3/2 field has the following 
form: 



(2.2) 



with 

A Q /3 = + M A )g af 3 



iA( la dp + 7/9 tt ) - -(3A 2 + 2A + 1)70^7/9 



M a (3A 2 + 3A + 1) 7q7/3 ], 



(2.3) 



where Ma is the mass of the spin | particle (A(1232) is the most familiar example), 
and A is an arbitrary parameter (A ^ — |). This Lagrangian is invariant under the 
transformation 
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A -2a 
1 + 4a ' 



(2.4) 



where a ^ — |, but is otherwise arbitrary. The Lagrangian 12.21 leads to the equations 
of motion 



AaB^ = 0. 



(2.5) 



The propagator for the massive spin 3/2 particle satisfies the following equation in 
momentum space 



K*p{k)G p s {k) = g aS , 



(2.6) 



where g a s is the metric tensor, and 



1 



Aa/3 = -[(-7 • k + M A )g a p - A{ la k !] + 1 pk a )-^A 2 + 2A + l) lal -k lp 



- M A (3A 2 + 3A + i) 7a7/3 ]. 



(2.7) 



Solving for G, one gets 
7 • k + M A 



G a0 {k) 



k 2 -M\ 
1 A + l 
3Ml 2A + 1 



9aP 3 7*7/3 



1 2 



3M 



lakf. 



A 



2A+1 



7ak(3 + 



1 A-l 

2 2A + 1 



7 • k 



\ 

AM A 
2A + 1 



(2.8) 



The physical properties of the free field are independent of the parameter A, so finally, 
with the parameter choice A — — 1, we obtain the Rarita-Schwinger [021 EH] propagator 
for the spin A particle, 



.^K + M A 



2k^k v 7^" -7" A;'*, 



(2.9) 



i A 3M A y 

This propagator, often found in the literature, is denoted by the double lines in the 
diagrams of Fig. 12.31 



We incorporate phenomenologically the effects of the non-zero width of the A by 
the replacement Ma — > Ma — zIa/2. This treatment of the finite width of the A is 
consistent with the Ward-Takahashi identities for the irpu vertex. This would not be 
true if Ta were introduced in the denominator of Eq. ()2.9j) only. 

Feynman rules 

In the vacuum, the nup vertex has the form 

-ilW = i^e^, (2.10) 

where g nup is a rrup coupling constant, and F n = 93 MeV is the pion coupling constant. 
We have also used a convenient short-hand notation 

e^ pQ = e^ a(3 p a Q p . (2.11) 

We choose Q as the incoming momentum of the tt, p as the outgoing momentum of 
the p, and q = Q — p as the outgoing momentum of the u, see Fig. 

P 



Figure 2.4: The assignment of momenta in the irup vertex. 

In the general case u and p are virtual particles. The value of g nujp can be obtained 
with the help of the Vector Meson Dominance Model from the anomalous 7177 decay. 
According to this model the electromagnetic interactions of hadrons are described by 
intermediate coupling to vector mesons, see Fig. 12. II The vector meson is subsequently 
converted into a virtual photon and finally the photon decays into the lepton pair, 
e + e~ . 

The vacuum 7r° — ► 77 decay 

One can calculate the g njl coupling constant from the width of the 7r° — > 77 decay in 
the vacuum. The transition amplitude of tt° — > 77 decay has the form 



(2.12) 



where the value of g wll = is model independent and results from Adler-Bell- 

Jackiw anomaly. The formula for the decay width is 

^^S 1 ^* 1 ^' (2 - 13) 

where, |p| = m n /2 is momentum of one of the photons in the rest frame of the pion, 
and the factor of 1/2 comes from the symmetrization over the indistinguishable bosons 
in the final state. Since the sum over all photon polarizations with momentum p gives 

a 

and from kinematics q 2 = m 2 , (Q — p) 2 = 0, we find 

£l-Mr = ^ 77 ^. (2-15) 



Consequently, the width has a form 

il 

TV 



2 m n 

1 7T77 = #7T7777Z> (2-16) 



where g^, applying the VDM (see Fig. 12. 1|) . can be expressed in terms of the g WLLip 
coupling constant as (see Fig. 12. 5|) 



e 



2 



^7T77 Zirup p2 n 2 n 2 

Finally, the width is 

2 S 

e m 



fi f 7T77 9nu>p 7-10 _o o ■ 

(2.17) 



^77 — g^ajp ni„2„2 raL- ( 2 -18) 



Using the physical parameters which we collect in Eq. 12.231 as well as T njl = 7.6eV, 
m n = 135MeV, e 2 = 47ra (where a = y|=), we obtain the value of coupling constant 
9-KLup by comparison of the 7r° — > 77 and 7r° — > pc<j diagrams, Fig. 12.51 It gives 
9-kujp = — ^ £ f ± fl , 7r77 = —1-36 (in our convention the sign of the coupling is negative). 



Meson-NN and meson-NA vertices 

The meson-NN vertices needed for our calculation have the standard form 



;9p /fj, _ lK p „pv^ \Jb 
2 2 771 j\r 



V b NN = " szr-^ft/K, (2.19) 



-iV*a NN = ^Ql 5 T a , 



2F n 
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Figure 2.5: Comparison of two diagrams, 7r — > 77 and 7r° — > up. Black dot indicates 

|r% coupling constants. The smaller black dot is equal to 



and 



„2 



^ or according to VDM. 

9u 9p ' & 



where p is the outgoing momentum of the virtual p, Q is the incoming four-momentum 
of the 7r, and a, b are the isospin indices of the tt, and p respectively, g u and g p are 
coupling constants for the u and p meson, k p is tensor coupling of the p meson, and 
gA denotes the axial coupling constant. 

For the interactions involving the A we use the same choice as in Ref. [60 : 



— iVNA a n a — 9NAnQ a T a , 
-iV NAap » = ig N A P ih 5 9 a "-l^ h P a )T\ 
-iV A -A^ = -^(7Y/-7V"-7Y ,1 + 7Y^ 



(2.20) 



where T a are the isospin \ 



transition matrices [63 . They are defined through 



the Clebsch-Gordan coefficients as follows: (|,/3|T M ||,i3) = (^lizp^l^-Iz) , with i 3 
and I3 denoting the isospin of the nucleon and A, respectively. In the Cartesian basis 
the explicit form reads 
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where the columns are labeled by 23 

3 



2' 2' 



(2.21) 

left to right, and the rows by I3 = 



3 1 _i 

2' 2' 2' 



.,, top to bottom. The following useful relation holds: 



rpa] rpb 



-6- 



ab 



-e abc r r . 



(2.22) 



3 3 

The values of the NA coupling constants follow from the non-relativistic reduction 
of the vertices and comparison to the non-relativistic values [B3| . The AAu; coupling 



incorporates the principle of universal coupling. It is obtained by replacing fc M — ► 
+ oj^ in the Rarita-Schwinger Lagrangian, Eq. 12.71 

Note that the pAA vertex does not appear in this study. Presence of pAA coupling 
in diagram 12 .HI implies that the uNA coupling should also appear. This is impossible 
from the Wigner-Eckart theorem, since u is an isoscalar particle. 

The choice of physical parameters is taken from other works [HI] and is as follows: 

g u = 10.4, g p = 5.2, ^ = 1.26, 

2.12 2.12v^2 
gNAn = , gNAp = • (2.23) 

In the results presented below we have chosen the tensor coupling of the p meson, 
k p = 3.7. The principle of universal coupling and the field current identity imply that 
the isovector part of the anomalous magnetic moment of the nucleon, Ky, is equal 
to the tensor coupling of the p meson, i.e. k p = Ky = 3.7, [65j, Ky is the isovector 
anomalous magnetic moment of the nucleon. However, the empirical value of k p / Ky 
extracted by Hohler and Pietarinen [HE] using a dispersion relation analysis of the irN 
scattering data is 1.78. Qualitatively, in our analysis similar conclusions follow when 
k p = 6 is used. 

We have to admit here that the relativistic form of couplings to the A, as well as 
the values of the coupling constants, are the subjects of an on-going discussion and 
research EZ1 ESI EHj ■ We are interested in estimating the size of the medium effect 
on the irup vertex and we do not search for precise and accurate predictions of tensor 
couplings. 

In order to simplify the approach we do not include any form factors in the vertices. 
In general, a vertex may depend on outgoing meson momentum squared, p 2 , q 2 , 
Fig. 12.41 and the Lorentz scalars, p • k and q • k, which enter the so-called sidewise 
form factors. Those form factors are not included in our calculations, because they 
lead to fundamental problems with gauge symmetry conservation. In particular the 
current conservation identities, Eq. 12.301 are violated. 

Tensor structure of the nup vertex in nuclear medium 

In our calculation we evaluate the diagrams of Fig. 12.31 and obtain the in-medium 

-Kujp vertex function 



where 



A'» u = A 1 e' lupQ + A 2 e» vkQ + A^ vpk + {Mtf + A^Q V + A & k v )e» kpQ 

+ (A 7 p» + A^ + A 9 k»)s vkpQ . (2.25) 

Thus in medium we have 9 possible tensor structures, while in the vacuum there was 
only one, e M ^ p< ^, Eg 12.101 All coefficients Aj are scalar functions of the four- vectors Q, 
p and k, with k° = E^. The occupation function is made explicitly Lorentz- invariant 



when we write |k| = a/ (k ■ u) 2 — m 2 N , where u is the four- velocity of the medium. 

The term with k 11 , upon the evaluation of the integral can be in general propor- 
tional to the three Lorentz vectors present in the problem, p^, Q^, and where w M 
is the four- velocity of the medium. 

/ WfE^ Am( ^ f ~ lkl)5{k ° ~ Ek) = Appfl + AgQfl + AuU ^ (2 ' 26) 

where A p , A g and A u are scalar functions of p 2 , Q 2 , p ■ Q, p ■ u, q ■ u, and fcp. 
Contracting Eq. 12.261 with Q^, p^ and u M we obtain a linear algebraic equations for 
A p , A q and A u which may be solved. To make this problem simpler we use the leading 
density approximation. This is equivalent to putting the three-vector k to zero in the 
functions A 2 — A g in Eq. 12.241 in the rest frame of the medium. Then we replace the 
loop momentum integration by 

/ w^ e(kf " miko ~ Ek) ~> h pB > (2 - 27) 

where p B is the baryon density. Now, with k = 0, k° = mjv, the contraction of 
Eq. 12.261 with g M , p^, and gives the set of equations 

1 ~p B m N p°A(k = 0) = A p p 2 + A q Q-p + A u p°, 

1 ~ PB m N Q A(k = 0) = App-q + A q Q 2 + A u q°, (2.28) 

— p B m N u°A(k = 0) = A p p ■ u + A q Q ■ u + A u u°, 

where u° = 1 in the rest frame of the medium. 

Since in the general case the vectors p, Q, and u are linearly- independent, the 
solution of Eqs. 12.281 is A p = A q = 0, A u = -^p B m N A()<L = 0). Thus, only the 
term proportional to m m in Eq. (|2.26|) is present in the leading-density approximation. 
Summarizing, above calculations, similarly as in [60], are equivalent to replacing k^ 
by tunu' 1 , what is valid for every occurrence of k in Eq. 12.241 For the terms in this 
equation involving second powers of k, the same prescription holds, as may be simply 
verified by an analysis analogous to the one presented above. 



Finally, we obtain the most general tensor structure of the form 



A*" 



= A^p® + A 2 e fluuQ + A 3 e^ pu + (A 4 p u + A 5 Q U + A 6 u u )e' 
+ {A 7 jf + AsQ 1 * + A 9 u^)e uupQ . 



fiupQ 



(2.29) 



This structure, restricted by the Lorentz invariance and parity, is more general 
than in the vacuum case Eq. 12.101 due to the availability of the four- velocity of the 
medium, u. The result of any dynamical calculation will assume the form Eq 12.291 
The vertex A^ u satisfies the current conservation identities: 



which serves as a useful check of the algebra. 

For the numerical analysis we have developed a program written in 'Mathematica'. 
In order to calculate traces from all diagrams of Fig. 12 .3| we have to apply the Feyn- 
mann rules, use Lorentz index contractions, and compute traces of the Dirac matrices 
(up to as many as 14). It is rather complicated from the technical point of view, 
thus we use FeynCalc [SI]- This is a package for algebraic calculation in elementary 
particle physics, which is very useful for many problems involving baryon loops with 
density-dependent nucleon propagators. All our calculations are covariant. 

2.2 Decays of particles at rest with respect to the 
medium 

We begin the presentation of our results for the case where the decaying particle is at 
rest with respect to the nuclear matter, q = 0. 

The in-medium 7r° — ► 77 decay 

First, numerical results are discussed for the process n° — > 77, because of its funda- 
mental nature. This is the decay into two real photons, where the pion is at rest with 
respect to the medium. In this case the four-vectors Q and u are parallel, and out of 
the 9 structures in Eq. 12.291 only A± and A 3 survive: 



Thus in the rest frame of the medium, where Q = (m n , 0, 0, 0), it is found that 



Q = 0, and p^A"" = 0, 



(2.30) 



A" u = A x e^ vQ + A z e^ pu . 



(2.31) 



AiQoe^ + A 3 u e^ p0 

(m n Ai + A 3 )e^ p0 = (A, + A 3 /m T )s^ pQ . 



(2.32) 



The following notation is introduced, which combines the vacuum and medium pieces: 



i e 2 



A^ = — (g^ + PbB) e^ pQ . (2.33) 

All the medium effects are collected in the constant B. The results of the medium- 
modified vertex at the nuclear-saturation density (p# = p = 0.17 fm -3 ) are presented 
relative to the vacuum value g vac = g^p^- In this analysis the medium effects are linear 
in the baryon density. We define 

g c s = gwpu + PbB, (2.34) 

with B denoting the medium contribution. Separate contributions to B coming from 
diagrams (a), (b) and (c) of Fig. 12.31 for the process tt° — > 77 are as follows: 



9a9p9uj( k p + !) 
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+5M A m^ + (4mjv + 3M A )m>^ - 16M A mjv - 11M A ), 

where Ma is understood to carry the width £T A - With parameters ()2.23|) we find 
numerically for the formal case r A = 0, 

B [a) = 97GeV~ 3 , B (b) = -3.1GeV~ 3 , B [c) = -102GeV~ 3 . (2.36) 

For the physical vacuum value of the A width, r A = 0.12 GeV, we find 

B {a) = 97 GeV" 3 , B {b) = -(2.6 + 1.5i) GeV" 3 , B {c) = -(87 + 42i) GeV" 3 , 

(2.37) 

and 

\g c s/g v , c \ 2 = |0.99 + 0.04z| 2 = .99, (2.38) 

whereas for nucleons only (diagram (a)) we would have 

IsS/Svad 2 = 0-91 2 = 0.82. (2.39) 

The imaginary values come from the term. From Eq. 12.381 and 12.391 we see that 
the effects of the A act in the opposite direction than the nucleons, resulting in an 
almost complete cancellation between diagrams (a) and (a+b+c) of Fig. 12.31 



The 7T° — > 77* decay 

The Dalitz decay of ir° into a photon and a lepton pair proceeds through the decay 
into a real and a virtual photon, and subsequently the decay of the virtual photon into 
the lepton pair, n° — > 77* — > r ye + e~ , see Fig. 12.61 The contributions to the constant 




Figure 2.6: The 7r° — > 77* — > 7e + e decay. Wavy lines indicates the p, uo or 7, dashed 
lines the pions, the hatched blob denotes the medium-modified vertex, the black dot 
is equal to or — - according to VDM. 

£? from all diagrams of Fig. 12 .HI for the processes n° — ► 7(^)7* and 7r° — > 7( p )7*, are 
listed in Appendix A. 

The virtual photon can be either isovector (p-type) or isoscalar (cu-type); thus 7^) 
denotes a photon which couples as an isoscalar and 7( p ) a photon which couples as an 
isovector. We denote the virtuality of 7* as M 7 * and investigate the dependence of g e & 
on M 7 *. The results both for the process it — > 7(^)7* and n° — > 7( p )7*, are displayed 
in Fig. 12.71 All calculations are done at the saturation density Pb = Po and with the 
vacuum value of the A width, Ta = 120 MeV. In the two upper plots of Fig. 12.71 
the solid lines show the real and imaginary parts of the full calculations, which means 
that we take into account all diagrams of Fig. 12. HI (with nucleons and delta isobars), 
whereas the dashed lines show the calculations with nucleons only which, of course, 
are real. The A resonance has imaginary contribution for ImT^, whereas for nucleons 
the result is always real. 

In the two bottom plots of the same figure, the solid and dashed curves indicate full 
calculations with deltas and nucleons, respectively. The results displayed in Fig. 12.71 
indicate that the change of My* in the allowable kinematic range from to has 
very little effect on the ratio of g e s/g vac and \g e s / g VSLC \ 2 ■ A somewhat different behavior 
for the isoscalar and isovector photons at M 7 * approaching is noticed. We observe 
that curves corresponding to the calculations with the A isobar shift the curve slightly 
up for 7T° — > 7(<j)7 decay, whereas for isovector photons the curves are pulled down, 
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Figure 2.7: The quantity g e s/gva. c and |g e ff/<?vac| 2 for the processes 7r° — > 7(^)7* and 
7t° — > 7( p )7*, plotted as a function of the virtual mass of 7*, M 7 *. The solid lines show 
the results of the full calculation with all diagrams (a+b+c) of Fig. 12. HI included, 
while the dashed lines correspond to the calculations with nucleon diagram (a) only. 
The upper plots include results of real and imaginary parts. All plots are done for 
the case ps = Po and Ta = 120 MeV. 



when approaching value. On the other hand for both of those decays the shape 
of curves without the A isobar is similar. By comparing the curves corresponding to 
the calculation with the A (diagrams (a+b+c) of Fig. 12. 3|) and with nucleons only 
(diagram (a)), it is visible that they act in opposite directions in diagrams (a) and 
(a+b+c), and it holds for the whole kinematically-available range. 

In Fig. 12.81 we show the result of the calculation with fixed r^., but with 
and Ma scaled down to 70% of their vacuum value. Such a reduced baryon mass is 
typical at the nuclear saturation density and is predicted by the Brown-Rho scaling 
hypothesis, see Chapter 1. The reduced mass of A is expected to behave similarly 
to the nucleon case. We notice that g e s is practically unchanged, and Fig. 12.71 and 
Fig. 12.81 show a similar behavior. 
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Figure 2.8: Same as Fig. 12. 7\ but for the values of and Ma reduced to 70% of 
their vacuum values. 

The uj — ► 7r°7* and p a — > n a Y decays 

Next, we consider the processes uj — > 7r°7* and p a — > 7r a 7*, where u and p are on the 
mass shell. In these cases the following vacuum values of the coupling constants are 
used, g^jm = —1.13 and g p ^-K = —0.76. They are calculated from the experimental 
partial decay widths 



Ta;7r7 = 717keV and r p7r7 = 79keV. 



(2.40) 



Let us concentrate on the first decay. The decay width for the process uj 
the form 



ujTvy 



4nM 2 p 2 p 



32vr 2 3m 2 ' 



e (l) £ (2) 



where 
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and with the momentum of the p 

Putting together all the needed parameters from Eq. I2.23l we obtain the physical value 
of the coupling constant, g^ p = —1.13. In a similar way we obtain g pnuj . Note that 
these coupling constants differ from each other, as well as from g nujp = — ^Jr- 5^-77 = 
— 1.36 used earlier for the 7r° decay. These differences may be attributed to form-factor 
effects: the virtuality of the vector mesons is changed from in the 7r° — > 77 decay to 
the on-shell value in the vector meson decay. For the case q = the contributions to 
B from Eq. 12.331 corresponding to the diagram (a,b,c) from Fig. 12 .3[ can be written 
as 

The formulas for Nu\ and are quite cumbersome, which reflects the presence of 
many terms in the Rarita-Schwinger propagator. We give them in Appendix B. 

There is an interesting phenomenon in the decays of vector mesons, related to the 
analyticity of the amplitudes of Fig. I2.3f b) and (c) in the virtuality of 7*. From 
the denominators of diagrams (b) and (c) of Fig. 12. 3[ see Appendix B for the case 
Ta = 0, we can see that the amplitudes from Eq. 18.31 have a pole at the value of M 7 * 
equal to 



Mq = rn v m? N + mlm N + m 2 v m N - M\m v + mlm v _ Q MQeV ^ ^ 

where m v = m p or m w . The above analytic structure is manifest in the numerical 
calculation presented below. For Ta 7^ the pole at 12.451 moves away from the real 
axis. Thus, analyticity is important and immediately leads to large changes of the 
vertex function near the poles. This formal feature of the presence of the pole can 
be clearly seen in the plots of Fig. 12.91 for the case of Ta = (dashed lines). For the 
values of Ta = 60 MeV (dotted line), Ta = 120 MeV, the physical value (solid line) 
and Ta = 240 MeV (dot-dashed line) the pole is washed-out, but its traces are still 
visible, with the curves reaching maxima around M = 0.4 GeV. The upper value 
of the virtual mass is M 7 * = m w — = 0.64 GeV. The differences between all 
the curves for various values of Ta in the range of M 7 * between 0.2 and 0.6 GeV 
shows that the results are sensitive to the assumed value for Ta- We can observe that 
at low values of M 7 * the effective coupling constant remains practically unchanged. 
However, at higher values of M 7 », above 0.2 GeV, the value of g e s is significantly 
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Figure 2.9: Top: \g e s/ '5w| 2 for the decay uj — > 7*7r° at the nuclear saturation density, 
plotted as a function of the virtuality of the photon, M 7 * . The dashed, dotted, solid 
and dot-dash lines correspond to the cases, Ta = 0, Ta = 60 MeV, Ta = 120 MeV, 
Ta = 240 MeV, respectively. Bottom: the same for the decay p a — > r j*Ti a . 



larger than in the vacuum. For the square of the coupling we find an enhancement 
by a factor of ~ 2 for the uj decay, and by a factor as large as ~ 5 for the p decay. 
All these numbers are quoted at the saturation density pB = Po- 

In Fig. 12.101 we display the results for different values of the baryon density, ps- 
For simplicity we denote the saturation density as 



Pb 

y = — • 

Po 



(2.46) 



We have done the calculations for both these decays for n p = 3.7 (dashed lines) and 
k p = 6 (solid lines). We can see that the effective coupling constant is more enhanced 
with growing density ratio y. 

Our method can be trusted numerically only at sufficiently low values of the baryon 
density, such that the leading-density approximation holds. For the masses or con- 
densates mentioned in Chapter 1 the nuclear saturation density its low enough to use 
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Figure 2.10: Top: \g e &/ '^vacl 2 for the decay u> — > 7*7r° at the nuclear saturation density, 
plotted as a function of the virtuality of the photon, M 7 « . The solid lines correspond 
to the case with n p = 6, dashed lines for n p = 3.7, y = Bottom: the same for the 
pa ^a^* decay. 



the low density approximation. On the other hand, we can see from Fig. 12.101 that 
the effects are large already at the saturation density. 

Therefore, numerical results at high densities obtained from our consideration are 
strongly indicative of possible large effects in the Dalitz decays of vector mesons in 
medium. In the uj — > 7r°7* and p a — > 7r a 7* decays, we compare the calculation with 
and without the A isobar, see Fig. 12.111 

In the figure we plot the ratio <7 e ff/<7vac and \g e s / g V a,c\ 2 as a function of the virtual 
photon mass, M 7 * at the saturation density y — 1 and for = 120 MeV, same as we 
have done in the previous case of the ir 77* decay. In two upper plots we show the 
imaginary and real (solid line) part of the calculations including nucleons and delta 
isobars, and the calculation with nucleons only (dashed lines). In Fig. 12.111 we can 
see that the real and imaginary functions of the considered processes, to — > 7*7r° and 
p a — > 7*7r a , are changed in comparison to Fig. 12.71 The features of the dashed lines, 
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Figure 2.11: The ratio of the effective coupling constant to its vacuum value at the 
nuclear saturation density for the decays uo — > 7*7r° and p a — > 7*vr a , plotted as a func- 
tion of the virtuality of the photon, M 7 * . The solid and the dashed lines correspond 
to the calculations with all diagrams of Fig. 12. 3\ and calculations with diagrams in- 
cluding only nucleons, respectively. Top plots show the real and imaginary parts for 
all considered cases. We use Ta = 120 MeV. 



corresponding to the calculations with nucleons only, remain practically unchanged. 
Both lower plots indicate that for the case with the A, the change of M 7 * in the 
kinematic range above 0.2 GeV has a visible effect on the effective coupling constant. 
We have already mentioned this when discussing Fig. 12.101 

Finally, we complete our discussion of the irup vertex at q = with Fig. 12.121 
which shows the result of the calculation for the values m N and Ma reduced to 70%. In 
this case there are significant differences compared to the unreduced mass, Fig. 12.111 
We can see that the effective coupling constant g e ff is enhanced by about 1.5 for the 
uj — > 7r°7* decay and by about 6 for p a — > n a ^* decay, and shifted to higher values of 
M r , above 0.4 GeV. 
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2.3 Decay of particles moving with respect to the 
medium 

So far we have presented our results for the nujp vertex for the decays of particles at 
rest with respect to the medium and we have pointed out large medium effects al- 
ready at the nuclear saturation density. Now, we will introduce the general kinematic 
case of finite three-momentum of the decaying particle with respect to the medium, 
and calculate the widths and spectral functions in the transverse and longitudinal 
polarization channels. 

In view of the substantial increase of the coupling constant for the uo and p decays 
when these particles are at rest with respect to the medium, we expect that the effect 
survives for particles moving with respect to the medium. We analyze the uo — > 7r°7* 
and p a —* 7r a 7* processes in the case where the uo and p meson move with a non-zero 
momentum with respect to the medium, denoted by q. The uo and p particles can 
have transverse or longitudinal polarization, defined by quantizing the spin along the 
direction of q. In medium the properties of these particles are different for these two 
polarizations; in particular, their widths T L and T T are different. 



Kinematics 

We consider the kinematics for the decay uj — > 7r°7* in the rest frame of the medium 
(not in the rest frame of the decaying particle as it is usually done in calculations in 
the vacuum). Thus the four- velocity of the medium is u = (1,0, 0, 0), and 

q 2 =ml, q-u = q = y/m* + |q| 2 , p ■ u = p , 
V ■ q = Poqo - cosa|p||q|, p 2 = p 2 - |p| 2 = M 2 *, (2.47) 

where q and p are the momenta of the u and 7* , respectively. Using the energy 
conservation: 



E w = E r + E n , (2.48) 



where 



E n = y/ml + jpp + |q| 2 - 2cosa|p||q|, (2.49) 



we find in general two solutions for the momentum |p|: 

cosa(m 2 r + M 2 » + 
2(M 7 2 * +q 2 sin 2 a) 



(i >2) ■ _ |q| cos a(m 2 +M 2 , + m 2 ) 



q Jmi - 2m 2 (M 2 » + m 2 ) + M^, + — M 2 « (4 sin 2 a|q| 2 + 2m 2 . 



p^> = y M 2 * + (pWy, (2.50) 

We consider two cases which are shown in Fig. 12.131 

Depending on the value of |q| we may have a case where only one branch for |p| is 
present, or where two branches appear. When |q| is smaller than some critical value 
<?crit we have the situation of Fig. 12. 13( A) . The values of |p| belong to one branch, 
and a takes all values between — 7r and tt. When |q| is larger than g cr ; t then we have 
the situation depicted in Fig. 12.13( B). For each value of a in the allowable range we 
have two branches. The angle a has a maximum value a* when |p x | = |p 2 |. Solving 
the equation corresponding to this condition we find that 



a = arcsm 



iryV^ " 2(M '* + ml)ml + (M '* " ml)2 ) ■ (2 ' 51) 



,2|q|M 7 » 

The transition between the behavior of Fig. 12.131 A) and B) occurs for q = q cr - lt . The 
value of g cr it follows from the condition |p 2 | =0 for a = as is clear from Fig. 12.131 
This condition yields, from Eq. 12.501 



mi - 2(M 2 , + m 2 )m 2 + (M 2 * - m 2 ) 2 

= fet. (2.52) 



2M 7 « 




Figure 2.13: The uj — > 7r°7* decay in the case where the uj meson moves with a non- 
zero momentum |q| with respect to the medium. Figure A): the situation when |q| 
is smaller than g crit , and the values of |p| belong to the first branch of Eq. 12.501 The 
angle a (between momentum |p| and |q|) takes values from -7T to tt. Figure B): for 
each value of a in the allowable range there are two branches. The second branch 
appears only for |q| above the critical value. 



Transverse and longitudinal polarizations 



Now we are equipped with appropriate kinematic expressions to consider the width 
of the transversely and longitudinally polarized u meson due to the decay uj — > 7r°7*. 
The expression for the width has the form 



, = -T- 

° 7 * n„ ^ 2o f 



d 3 p 



d 3 p' 



T \M ss/ \ 2 (2n) 4 5^(q- P ~p'), (2.53) 



J2q j {2ixf2p J (2n^2p 

where s is the polarization of the uj, s' is the polarization of the 7*, n s is the number 
of spin states of the uj meson, and |.M SS '| 2 is the transition amplitude. We perform 
the phase-space integral with the kinematics discussed above and obtain 
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where Y2b * s ^ ne sum over ^ ne one or ^ wo branches, and expressions 
(1 2) ^(^o — V M 2 + r 2 — \Jm 2 + r 2 — 2r|q| cos a + q 2 ) 



ch- 



are the factors from the Dirac delta function Si^E^ — E p — E n ). 



(2.55) 



= |p(l,2)| 



The uj meson has transverse polarization described by the polarization vectors 
£^_n , and £^ n . Using explicit forms of these polarization vectors one finds the following 
formulas [HI EE]: 



' £ (+f(+) - £ (~) £ (~) 



q ■ q — (q ■ u) 1 



(2.56) 

q ^ qU . ,w + (t-q-uu»){q» -q-uu») _ 
~ £ (o) £ (o) — \-uu H — = L [q , 

The tensors and are defined with such signs as to form projection operators, 
i.e. , T^ u T u a = T^, L^L? = L^ a , and T^L? = 0. Furthermore, we have T^q v = 
and L^ v q v = 0, which reflects current conservation, as well as T^ v u v = 0. 
By summing over all polarizations, we find the usual expression 

- E e WM e (-> = + ^ = <T - — = Q" v (q)- (2-57) 

s=±,0 ^ ^ 

Now we return to an expression for the width, Eq. 12.531 We write the transition 
amplitude 

ss' ss' 

fig 

s 

where £( s /w is a polarization vector for 7* and E( s \ v is the polarization vector for 
the uj meson. Summing over s = 0, and s = ± we obtain, respectively, 

\M L \ 2 = LftQftfaA^ 

\M T \ 2 = l^Q^A^A*^,. (2.59) 
Combining Eq. 12.541 and 12.591 we may write the final expression 
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Results for transverse and longitudinal widths 

Our numerical results for the uj — > 7r°7* decay are shown in Fig. 12.141 We present the 
dependence of widths for the transverse and longitudinal polarization on the virtuality 




0.1 0.2 0.3 0.4 0.5 0.6 



M r [GeV] 

Figure 2.14: The longitudinal (top) and transverse (bottom) widths for the uo — > 7r°7* 
decay, plotted as a function of the virtuality of the photon, M 7 . . The dashed and 
solid lines show the vacuum and medium cases for ps = po, respectively. 

of the photon, M 7 *. We have used here the constant Ta = 120 MeV. Calculations are 
done for different values of |q|. From Fig. I2.14l we can conclude that the transverse and 
longitudinal widths decrease with momentum |q|, and the medium effect is weakened 
with increasing momentum |q|. However, we can observe that the effect remains 
substantial for |q| up to ~ 200 MeV, the relevant value for temperatures typical 
in heavy-ion collisions. As it is commonly known in a fireball formed in relativistic 
heavy- ion collisions, the momenta are lower than a typical temperature ~ 150 MeV. 

We present analogous results for the p a — > 7r a 7* decay in Fig. 12.151 We notice 
again that while |q| increases from to 600 MeV, the medium effect is weakened, 
and both the transverse and longitudinal widths decrease. 

We see that the medium effect is large for the typical momenta up to ~ 200 MeV. 
For higher momenta this effect is lower; however it still remains visible compared to 
the vacuum result, see the dashed line in Fig. 12.151 
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Figure 2.15: The longitudinal (top) and transverse (bottom) width for the p a — > 7r a 7* 
decay, plotted as a function of the virtuality of the photon, M 7 *. The dashed and 
solid lines show the vacuum and the medium cases, respectively. 



Chapter 3 
Summary of Part I 



The main conclusion of our investigations in this part of the thesis is that the medium 
effects on the irup coupling in the discussed Dalitz processes uo — > 7r°7*, p a — > 7r a 7* 
are large. Moreover, these medium effects come dominantly from the processes where 
the A isobar is excited in the intermediate state. On the other hand, for the tt° — * 77* 
the effects of the A isobar cancel almost exactly the effects of the nucleon particle- 
hole excitations, such that the medium effect is small. The Dalitz yields from the 7r° 
decays are not altered by the medium for two reasons: first, virtually all pions decay 
outside of the fireball due to their long lifetime; second, the width for this decay is 
practically unaltered by the medium. 

At the nuclear saturation density, for the Dalitz decays from the p and u mesons, 
the value of the effective coupling constant is enhanced compared to the vacuum value. 
The increased coupling constant greater by about a factor of 2 for the uj meson and 
by about a factor of 5 for the p, leads directly to large widths of the uj and p mesons 
in medium. 

Next, we have analyzed the case where the decaying particle (p or a;) moves with 
respect to the medium. In this case we have to look separately at the longitudinal 
and transverse polarizations. The properties of mesons are different for these polariza- 
tions. We investigated the dependence of the widths for transverse and longitudinal 
polarization on the invariant mass of dileptons. This leads to the conclusion that the 
medium effects decrease with growing magnitude of three-momentum, but remain 
significant for typical momenta in heavy-ion collisions, q < 200 MeV. 

The enhancement of the coupling constant directly affects the calculations of the 
dilepton production in relativistic heavy-ion collisions. An increased value of width 
results in an increased dilepton yield. In the next part we will apply our model to 
evaluate the dilepton production from the p and uj decays in heavy-ion collisions. 
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Part II 
Dilepton production rate 
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Chapter 4 



Electromagnetic signals from hot 
and dense matter 

In this part of the thesis we use directly the results of Part I. We focus on the pro- 
duction of dileptons in hadronic decays, which is an area of great interest from a 
theoretical point of view. We will investigate the theoretical yields from the Dalitz 
decays of vector mesons in the region of 0.2 — 0.6 GeV, where many existing cal- 
culations have problems to explain the CERES and HELIOS experiments jU EJ 
E21 EHl EH EE1 EHl EE! • O ur purpose is to see to what extent the medium modifica- 
tion of the nujp vertex influences the dilepton production, and may be significant in 
understanding the experimental data in the region of low invariant masses. 

As we mentioned in the Introduction, electron pairs are produced in all stages of 
ultrarelativistic heavy-ion collisions. Because of the small electromagnetic coupling 
constant, they carry undistorted information about the densest and hottest phases of 
the heavy- ion collisions where chiral symmetry restoration and/or quark gluon plasma 
(QGP) formation are expected. 

The first pairs of leptons are created during the stage when two nuclei pass through 
each other, while the last pairs are formed when the hadrons are decoupled and 
move freely to the detectors. After decoupling of hadronic matter the greater part of 
dileptons is produced by the electromagnetic decays of hadronic resonances. For pair 
masses below 140 MeV/c 2 the tt° Dalitz decay dominates the spectrum, whereas at 
higher masses the u, p, r], rj and decays should be most significant. 

The dilepton mass spectra are very important for the study of the in-medium 
properties of vector mesons, which we have discussed in detail in Part I. 

As we have pointed out in the Introduction, experimental measurements of dilep- 
ton mass spectra can be divided into three mass regions. The low mass region, below 
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1 GeV, is dominated by hadronic interactions and hadronic decays at freeze-out. It 
is particularly sensitive to in-medium modifications of the light hadrons. The inter- 
mediate mass region lies between 1 GeV and about 2.5 GeV, where the contribution 
from the thermalized QGP might be seen. We want to mention that the physical 
picture of QGP structure has changed in the last year. Now, many theoretical devel- 
opments [ZZIIZH1 suggest that the matter produced at RHIC, in the high temperature 
region T c < T < 4T C , is a strongly coupled quark gluon plasma (sQGP) contrary to 
the previous expectations which were based on weakly coupled quasiparticle gas. In 
order to understand what sQGP is, one may have a look at two other examples of the 
strongly coupled systems such as trapped ultra-cold atoms in a large scattering length 
regime or supersymmetric Yang-Mills gauge theories at strongly coupling limit (Con- 
formal Field Theory with a non-running coupling at finite temperature). Finally, the 
high mass region at and above mj/^ is important in connection with the suppressed 
J/ty production with respect to the background from the Drell-Yan process of the 
quark- ant iquark annihilation. 

So far, the experimental measurements of dilepton yields in heavy-ion collisions 
have mainly been carried out at the CERN SPS facility by three collaborations: the 
CERES (Cherenkov Ring Electron Spectrometer) collaboration has measured dielec- 
tron spectra in the low- mass region @JE3E|, the HELIOS-3 collaboration has spe- 
cialized in dimuon spectra up to the J/ty region 0, and the NA38/NA50 [TH } ITT ^ IT7] 
collaborations have been dedicated to study of dimuon spectra in the intermediate and 
high-mass region, emphasizing the J/^> suppression. Dilepton data have also been 
taken by the DLS collaboration at BEVALAC ED] at much lower bombarding 
energies of about 1 — 2 A GeV. 

The study of the low mass dileptons between 0.2 — 0.6GeV has recently become of 
considerable interest. In this region the dilepton continuum originates from the Dalitz 
decays of neutral mesons such as ir°, 77, rf — > e + e _ 7 and u — > e + e~7r°. The resonance 
peaks are due to direct decays, u, p, <fi — > e + e~ . In central heavy-ion collisions the low 
mass dilepton enhancement is observed by the CERES and HELIOS-3 collaborations. 
This enhancement has been studied in various approaches, ranging from hydrodynam- 
ical models to complicated transport models. The calculations based on the cocktail 
model, which are designed to describe the proton-induced reactions, work very well for 
light-heavy systems, but fail to explain the observed phenomenon in nucleus-nucleus 
collisions. Various medium effects, such as the dropping vector meson masses, as first 
proposed by Brown and Rho £Q, the modification of the rho meson spectral func- 
tions, or the enhanced production of 77 have been proposed to explain the observed 



enhancement. 

In the intermediate-mass region from about 1 GeV to about 2.5 GeV the excess 
of dileptons has also been observed by both the HELIOS-3 and NA38/NA50 collabo- 
rations. This enhancement is observed in central S+W, S+U, and Pb+Pb collisions 
in comparison to proton- induced reactions E3 El El- The intermediate- mass 
dilepton spectra are particularly useful in the search for QGP. 

The J/^f suppression in the high- mass region above 2.5 GeV is a very interesting 
problem. It was interpreted as a signal of QGP formation by Matsui and Satz already 
in 1986 |79j . Nowadays, there is an indication from experimental data that up to 
central S+Au collisions pre-resonance absorption in nuclear matter is sufficient to 
account for the observed J/^f suppression [HU]. However, data for central Pb+Pb 
collisions from NA50 collaboration show an additional strong J/ty suppression, which 
is interpreted gnal of color deconfinement [8*T) l%2~] . 

Below we discuss mainly the enhancement of the low-mass dileptons, because we 
have found, in Part I, that the considered nup coupling constant is sizeable enhanced 
in the medium just in the region between 0.2 — 0.6 GeV of the invariant mass of 
dileptons. Hence in this region the modifications of the vertex directly influence the 
dilepton production rate. 

4.1 Overview of experimental measurements 

The CERN SPS is the first machine used for systematic investigations of the dilep- 
ton production in ultrarelativistic hadron-nucleus and nucleus-nucleus collisions. In 
this section we review the CERN-SPS dilepton experiments CERES, HELIOS, and 
NA38/NA50. Since the dilepton spectra are usually compared to the expected con- 
tributions from hadron decays, the co-called hadronic cocktail, we begin this section 
by a short description of the hadronic cocktail model. 

Hadronic cocktail 

The hadronic cocktail is used to simulate the relative abundance of dielectrons pro- 
duced by hadron decays in proton-nucleus and nucleus-nucleus collisions at the in- 
variant mass range covered by the CERES acceptance (m ee < 2 GeV) . This invariant 
mass range is dominated by the decay of light scalar and vector mesons, such as 7r°, rj, 
rj', p°, u, and (p. The hadronic cocktail treats proton-nucleus and nucleus-nucleus col- 
lisions as a superposition of individual nucleon-nucleon collisions and gives a reference 
to compare with the experimental observed dilepton yield. As input for this simula- 



tion one needs to know differential cross sections, the widths of all decays including 
dileptons in the final state, and the decay kinematics for all involved particles. 

The proton-nucleus collisions are approximated by a superposition of nucleon- 
nucleon collisions, i.e. the proton-nucleus yield is obtained by rescaling the nucleon- 
nucleon yield with the mean charge-particle multiplicity of the colliding system. For 
Pb-Au collisions the relative cross sections (a/a n o, where o is the cross section for a 
particular hadron decay into the electron pairs) are taken from a thermal model, [H3J - 
In order to compare with experiments the hadronic cocktail contribution is divided by 
the total number of charged particles within the acceptance of the detector. The sum 
of all contributions from hadron decays is calculated with Monte-Carlo event genera- 
tor. The invariant mass spectrum is calculated including all known hadronic sources. 
The particle ratios of these sources are assumed to be independent of the collision 
system and to scale with the number of produced particles. Their p^-distribution are 
generated assuming m_i_-scaling, [S3] based on the pion p± spectrum from different 
experiments and fitted [85J to the formula: 

dN m_L 

— ~ P± exp[-— ], (4.1) 

where for SPS at energy 158 A GeV the inverse slope parameter is parameterized as 

T eff = 175 GeV + 0.115 m[GeV], (4.2) 

with meson mass indicated as m. The rapidity distribution is a fit to measured data. 
All Dalitz decays are treated according to the Krall-Wada [HE] expressions with the 
experimental transition form factors taken from Ref. [87J. The vector meson decays 
are generated using the expressions derived by Gounaris and Sakurai [88.. 

Charm production is not taken into account since it is negligible in the low-mass 
range. Finally the laboratory momenta of the electrons are constructed with the 
experimental resolution and acceptance. 

The differences between hadron cocktail and experiment indicate the existence of 
in-medium effects or the violation of the scaling behavior. 



4.1.1 CERES/NA45 experiments 

The only experiment which measures the low- mass dilepton pairs up to about 1.2 GeV 
is the CERES/NA45 experiment. This collaboration carries out systematic measure- 
ments of dilepton spectra in proton-induced reactions. Figure l4~T| taken from [89] . 
shows the data on dielectrons of proton-induced collisions at 450 GeV in comparison 
to the expected contributions of hadron decays, the hadronic cocktail For instance, 
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Figure 4.1: Dilepton invariant mass spectrum of p+Be and p+Au collisions at the 
incident proton energy of 450 GeV, taken from [HH]- The expected sources from 
hadron decay (solid lines) and the experimental points are in excellent agreement. 

the dashed curves in Fig 14. II give the expected dielectron spectra from hadron decays 
for p-Be and p-Au reactions. The total expected cocktail of hadronic sources is de- 
noted by the solid line. We observe that the dilepton yield is fitted by the cocktail 
model very well. The situation is different for dilepton spectra in central S+Au and 
Pb+Au collisions at 200 A GeV and 158 A GeV, respectively, see Fig. l4~2l from 
The dielectron yield significantly exceeds the expectations extrapolated from proton- 
induced reactions. The enhancement is most pronounced in the invariant mass region 
0.2 — 0.6 GeV. In both cases the difference between the experimental and the cock- 
tail model is about a factor of 6, whereas in p-induced reactions the spectrum from 
0.2 to 0.6 GeV reproduces the data. The CERES collaboration has also studied the 
low-mass dielectrons in Pb+Au collisions at the lower energy of 40 GeV, see Fig. 14.31 
Enhancement of the dilepton yield, related to the expected yield of hadronic sources, 
is even more pronounced than at the energy of 158 A GeV |5D] . 

The difficulty in explaining this experimental fact has become known as the dilep- 
ton puzzle and is an important theoretical or experimental problem. 

4.1.2 HELIOS experiments 

In the HELIOS experiments there are three collaborations. The HELIOS-1 collabo- 
ration jHI] was the pioneer in the measurement of the dilepton spectrum at CERN 
SPS. It was the only experiment that measured both the e + e~ and pair produc- 
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Figure 4.2: Dilepton invariant mass spectrum of S+Au and Pb+Au collisions at 
200 GeV and 158 GeV from CERES collaboration, 89J. The data is compared to 
expected cocktail of hadronic sources. A visible increase of dilepton production is 
observed in the region 0.2 — 0.6 GeV. 



tion in p+Be reactions at 450 GeV as well as the first in which the hadron cocktail 
contribution to the dilepton spectra was analyzed. The conclusion is that the mea- 
sured dielectron and dimuon spectra in p+Be collisions are described by the cocktail 
hadron decays, which was later confirmed by the CERES collaboration with much 
greater precision. The HELIOS-2 collaboration [OH] was the first to measure the pho- 
ton signals using the S and O beams at the CERN SPS energy. Very small direct 
photon enhancement (of about a few percent) was observed in these studies. This 
observation was later confirmed by the CERES and WA80 collaborations with better 
statistics. 

The HELIOS-3 experiment |3] is designed to study the dimuon production in p+W 
and S+W collisions at 200 GeV, and covers the low- and intermediate-mass region. 
In Fig. 14 A\ taken from [02], we see the results for p+W collisions (left) and S+W 
collisions (right). For p+W collisions the low mass dimuons are mainly from the 77 
Dalitz decay, while the omega meson Dalitz decay is important in the mass region 
0.2 — 0.6 GeV. The cocktail contribution is in agreement with the experimental points. 
For the S+W collisions the dimuon production is enhanced in both the low-mass and 
intermediate-mass regions, as compared to the proton-nucleus collisions. 

In the low dimuon invariant mass region the enhancement factor is smaller than 
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Figure 4.3: Dielectron invariant mass spectrum of Pb+Au collisions at 40 A GeV 
taken from [89]. The data is compared to the expected cocktail of hadronic sources 
(solid line). 



that found by the CERES collaboration in central collisions. One of the possible 
reasons for this difference is that the HELIOS-3 experiment covers more forward 
pseudorapidities, 3.7 < rj < 5.2, than the CERES experiment which is close to the 
midrapidity region, 2.1 < rj < 2.65. A significant excess in the dimuon yield in the 
intermediate-mass region is also reported, with an enhancement factor of about 2 [S|. 

4.1.3 NA38/NA50 experiments 

The NA38 and NA50 experiments specialize in the study of the J/^f suppression in 
heavy-ion collisions. These collaborations measured the dimuon spectra in a mass 
region from the threshold to about 8 GeV. The measurements are carried out at 
energies of 200 A GeV for the p+W and S+U reactions (NA38), and 160 A GeV 
for the Pb+Pb reactions (NA50). A situation similar to that described for other 
experiments occurs. 

In proton-induced reactions the data are described with a rather good agreement 
as compared to the expected sources from hadron decay. On the other hand, for S+U 
and Pb+Pb reactions significant enhancement for low- and intermediate-mass region 
is observed. This shift is more visible in Pb+Pb collisions in the mass region around 
2 GeV. The data from NA38 /NA50 experiments have no centrality selection, whereas 
the CERES collaboration measures central collisions. 
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Figure 4.4: Dilepton invariant mass spectrum of p+W (left) and S+W (right) colli- 
sions at 200 GeV from HELIOS-3 collaboration. The data is compared to expected 
cocktail of hadronic sources, jH2] • 



4.1.4 Dilepton spectra at BEVALAC energies 



It is worth mentioning that the measurements of the dilepton production at much 
lower bombarding energies, in the range 1 — 5 A GeV, have been performed by the 
DLS collaboration at LBL BEVALAC [T9~l 120j. In this experiment a different tem- 
perature and density regime is probed. The dielectron spectra were measured for 
p+Be collisions at 1, 2.1 and 4.9 A GeV, for Ca+Ca reactions at 1 and 2 A GeV and 
for Nb+Nb at 1.05 A GeV. In the low invariant mass region from 0.2 — 0.6 A GeV 
production of dielectrons is found to be enhanced compared to estimates based on 
the theoretically known dilepton sources. 

For ultra-relativistic reactions (SPS) the low-mass dilepton excess is explained by 
the reduction of the p-meson mass in a dense medium. The in-medium broadening 
of the p-meson is also expected to be sufficient to explain the dilepton yield at SPS 
energies. These explanations fail for the DLS data. This fact is known as the 'DLS 
puzzle'. The reason probably lies in the fact that the possible pQCD contributions 
and sufficient -k + -k~ annihilations are absent at this energies. 



4.1.5 Other experiments 

The new NA60 experiment studies the production of open charm and the production 
of dimuons instead of dielectrons in proton-nucleus and nucleus-nucleus collisions at 
CERN SPS. This experiment is a continuation of the NA38 and NA50 experiments 
but with significantly enhanced detector technology. The aim of the NA60 experiment 
is to examine various possible signatures of the transition from hadronic to deconfined 
partonic matter e.g., charmonium suppression, dimuons from thermal radiation, and 
modifications of vector meson properties. The NA60 experiment is continuing with 
115 In beams at least up to the year 2005, Ref. [THj . 

Dilepton spectra were also measured at KEK in p+A reactions at a beam energy 
of 12 A GeV [94j . Also here considerable enhancement over the expected yield from 
hadronic decays was observed below the p-meson peak. Explanations of the experi- 
mental spectrum within a dropping mass scenario or the broadening of vector meson 
are not proper in this they were not at the Bevalac energies. Similar experi- 

ments with much better statistics are planned at GSI SIS by the HADES collaboration 

Summarizing this section, in the low-mass region between 0.2-0.6 GeV of the 
invariant mass of the dilepton, experiments reported an unexpected enhancement of 
dileptons in nucleus-nucleus collisions with comparison to the extrapolation from more 
elementary proton-induced reactions. This has been a topic of great interest from a 
theoretical point of view, to which the next section is devoted. 

4.2 Overview of theoretical models 

In this section we would like to review various theoretical explanations that have been 
put forward to explain the experimental data discussed in the previous section. A lot 
of theoretical efforts has been devoted to understanding the experimental results in the 
low-mass region from the CERES and HELIOS experiments. Theorists have proposed 
possible explanations of the dilepton puzzle based on transport model calculations, 
many-body approaches, and the effects of dropping rho meson mass and/or increasing 
width. 

The main contributions to e + e~ pairs with mass below 300 MeV are: the direct 
leptonic decay of vector mesons such as p°, uj and <p, the pion-pion annihilation which 
proceeds through the p° meson, the kaon-antikaon annihilation that proceeds through 
the <p meson, and the Dalitz decay of tt°, n, uj and p. 

The CERES and HELIOS-3 data are in reasonable agreement with each other in 



proton-induced reactions. They are well explained by the Dalitz and direct vector 
meson decays. In nucleus-nucleus collisions the substantial enhancement of low-mass 
dileptons in the mass region 0.2 — 0.6 GeV remains unexplained. The difference 
between proton-induced and nucleus-nucleus reactions may lie in additional contribu- 
tions from pion-pion and kaon-antikaon annihilation to dilepton production in heavy- 
ion collisions. Authors of (2H E21 El E3] have included the contributions from the 
pion-pion annihilation but neglected any medium effects. 

The theoretical results are still in disagreement with data in the mass region of 
interest, which leads to the suggestion that medium modifications of vector mesons 
are needed [3U EHl EH] ■ It has been found that the models with the use of in- medium 
vector meson masses describe the CERES data much better than models with vacuum 
meson masses. The same situation occurs for the dimuon spectra from the central 
nucleus-nucleus collisions measured by the HELIOS-3 collaboration [§21 EH]- 

The rho meson decays dominantly into two pions. The medium modifications of 
the pion dispersion relation and the rho-nucleon scattering contribute to the spreading 
of the rho meson spectral function, which results in enhancement of the production 
of low-mass dileptons from in-medium rho meson decay. This broadening of the rho 
meson spectral function provides a better description for the CERES (HEl EE] and 
HELIOS-3 data. 

Many calculations and mechanisms have been proposed for the explanation of the 
low-mass dilepton enhancement. Kapusta et al. |HZ| and Hung et al. [SB] proposed 
that the in- medium reduction of the rj and/or rj' masses leads to an enhanced produc- 
tion of these mesons, which in turn increases the low-mass dilepton yield. The authors 
of 1001 HOP! studied the effects of an in-medium pion dispersion relation on dilepton 
production and they indeed found that the dilepton yield is enhanced near 0.4 GeV 
invariant mass. The contribution from the pion-rho annihilation has been studied in 
|101[ I102j . This contribution was too small to explain the observed enhancement. 
Hung and Shuryak [103 also emphasized the effect of dropping a\ meson mass, which 
is the chiral partner of the rho meson, on the very low-mass dileptons. On the other 
hand, Kluger et al. |104j studied dilepton production rates from the disoriented chiral 
condensate (DCC) state, and found that the dilepton yield with mass below 300 MeV 
from the DCC state can be significantly larger than that from the equilibrium state. 

In other words, the low-mass dilepton enhancement reported by the collabora- 
tions from the previous section required the introduction of pion-pion annihilation 
and medium modifications of the meson properties. The dropping rho meson mass 
scenario, as implied by the Brown-Rho scaling conjectures and QCD sum rule calcu- 



lations, provides a possible explanation of the observed enhancement. On the other 
hand, the medium modification of the rho meson spectral functions due to its inter- 
action with baryons and medium modification of coupling constants also provides a 
reasonable description of experimental data. 

Here we have discussed various theoretical approaches. Next in this thesis we 
follow the lines of Brown- Rho pQ and Rapp-Wambach [HEl EH] and analyze medium 
modification of vector mesons. 



Chapter 5 



Formalism of the dilepton 
production from vector mesons 

We are interested in the nup coupling constant which enters into the Dalitz decay of 
the uj and p meson. The Dalitz p decay is usually neglected in other analyses because 
of its absence in the Particle Data Book |l()5j . Nevertheless, it may occur and while 
it is less significant in the vacuum than the uj Dalitz decay, we will show that in the 
medium it becomes equally important. 

We point out that we work at zero temperature. The authors of Ref. have 
analyzed the finite temperature effects for the pmr coupling. Such effects turn out to 
be very small, therefore in our analysis they are also neglected. 

In our calculations we use the Vector Meson Dominance Model to estimate the 
dilepton yields from vector-meson decays. Apart from the Dalitz decays, an important 
source of additional dileptons in the region near twice the pion mass is the tttc anni- 
hilation. In this Chapter we first briefly review the VMD model. Next, we describe 
the structure of three-body Dalitz decays. 

In order to calculate the dilepton spectrum one has to assume the evolution model 
of the collision system. Therefore, we review the Rapp-Wambach 12ZI model of 
the hydrodynamic expansion of the fire cylinder, which we are going to apply in 
our calculations. In section 4.4 we discuss the kinematics for both direct and Dalitz 
decays with inclusion of effects of medium expansion in heavy-ion collisions. Thus we 
get the yield of dileptons produced during the expansion. In this analysis the Bose 
enhancement is neglected for brevity. Finally, we present our results and compare 
them to the CERES experimental data. 
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5.1 Vector Dominance Model (VDM) 

The dilepton production rate per unit time and volume from the annihilation of pions, 
Fig. 15. 1[ is given by the formula 



7T + \ 




Figure 5.1: Vector Dominance Model of dilepton production from pion annihilation. 
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The sum runs over spins of leptons, where p\ and p 2 are the four-momenta of the 
annihilating pions and q± and q 2 are the four-momenta of the outgoing leptons. The 
dP and dQ indicate the integration measures 

dPi = w^E~: dQi = (2tt)32v (5 - 2) 



where E Pi = \/pj + ml and E q . = y/qf + mf, i = 1,2. We use the following normal- 
ization for the lepton spinors, uu = vv = 2m;. In Eq. 15.11 /(p) is the pion distribution 
function of the form 



exp 1 



(5.3) 



with the pion chemical potential fi n , incorporated in several works. The positive pion 
potential increases the number of pions over the equilibrium value, as discussed later 
on in section 5.3. 

The matrix element Ai^+^-^i+i- is calculated in the framework of the VDM. 
According to this model the electromagnetic interactions of hadrons are described by 
the intermediate coupling of hadrons to vector mesons, Fig 15.11 The annihilation of 
the pions (7r + 7r~ — > can be understood as the formation of a vector meson by 

the annihilation of pions and, the sequential decay of the vector meson into a lepton 
pair. Therefore, pions regain an electromagnetic form factor. The matrix element for 



this process has the structure 



M 



K+n-^i+i- = 4ei/ V y —g n nv(Pi - P2) » (5-4) 

<T y Jvm^-q 2 - im v Y v 



where the term m v /(m v — q — imyTy) comes from the vector meson propagator 

g^v _ q^q 

Gy(q) = 2 ^ . (5.5) 

q z — m v — imyL y 

The quantities my and Ty are the mass and width of the vector meson V, the coupling 
of the vector meson to the virtual photon are described by the fy parameter, g n7T y 
is the effective coupling constant for the 7nrV interaction, and 1/q 2 comes from the 
photon propagator such as 

„u _ 

with q denoting the four- momentum of the virtual photon (q 2 = M 2 , the invariant 
mass squared of the lepton pair). 

The quantity L M is the leptonic current which may be expressed in terms of the 
Dirac bispinors, 

L M (<7i, sx; q 2 , s 2 ) = 2m^(gi, si)j^v(q 2 , s 2 ), (5.7) 
where Si and s 2 denote the spin variables of the outgoing leptons. The term 

2 

EG i^y / \ / \ 

T^ 2 — „2 im r 9«*v(j>i-P2h (5.8) 
v Jv my — q — imyl y 

from Eq. 15.41 describes the coupling of the photon to two incoming pions via the elec- 
tromagnetic form factor in the VDM including also the propagator of the intermediate 
vector meson. 



5.2 Dalitz decays 

The structure of the three-body Dalitz decays of the form A — > Bl + l~ is as follow: a 
particle decays A into B and a vector meson V which is subsequently converted into 
a virtual photon according to the vector dominance principle. Finally, the photon 
decays into the e + e~ pair, see Fig 12.11 This process is usually called the 'internal 
conversion' of a virtual photon into a lepton pair. For these processes, where either a 
vector meson (V) or pseudoscalar meson (P) decays via V — > Pl + l~ or P — > Vl + l~ , 
the transition amplitude has the invariant form 

Ma^biH- = L^ a ^p a q^f AB (M 2 ), (5.9) 



where L M is a leptonic current, p a the four-momentum of meson B, q the four- 
momentum of the virtual photon 7*, q 2 = M 2 the invariant mass squared of the lepton 
pair, f A s{M 2 ) is the electromagnetic form factor describing the A — > B transition, 
and e a ^ lfl is a totally antisymmetric tensor. We easily obtain the useful decomposition 

\M A ->m+i-\ 2 = \M A -*bA 2 x v^iI-Mt^+H 2 , (5.10) 
whereby we can see why this process is called the internal conversion. Here 

\M r ^ l+l -\ 2 = e 2 L^L»'y (5.11) 

and 

\M A -+Br\ 2 = M A n'Y> (5.12) 

where 

A, = e al3 ^p a q^ef AB (M 2 ). (5.13) 

Using the decomposition E3D for the differential decay distribution we can write 

dM 2 = "Tm^^*' (5 - 14) 
where for the width of the internal conversion of the photon we find 



r r ^ = 2 -^fH ] /l-^L(2m l + M% (5.15) 



and 

r A ^ Br = r A ^ B ,\f AB (M 2 )\ 2 V, " A l ""'7; L "^,„ v """ (5.16) 



r2 ,12 i m A ~ (™>b + M) 2 f' 2 [m 2 A - (m B - Mff/ 2 
(m\ — m 2 B ) 3 



is related to the decays where the photon is real. The normalized form factor for the 
A — > B'-f* transition has a form 

o o m 4 n + m 2 T 2 

\Sab{M )\ = - ' ' ' (5.17) 
{nip — M A y + mil A p 

It has the property \f AB {M 2 )\ 2 -> 1 for M -> 0, for increasing M the form factor 
also increases. After putting all Eqs I5.14j I5.15j 15.161 and 15.171 together we arrive at 
the differential decay distribution which we will use to calculate the mass spectrum 
of dileptons from Dalitz decays. For more details see 106 . 



5.3 Model of the hydrodynamic expansion of the 
fire-cylinder 



The description of medium effects on the dilepton production in heavy-ion collisions 
requires modeling of the time evolution of the collision system. The theoretical models 
can be divided into: transport models [92 ^ I107j . hydro dynamical approaches 1IU3| I108( 
I102| . and thermal fireball models |36| l26j. 

In this work our calculations are based on the fire cylinder hydrodynamic expan- 
sion model of Ref. |2T^I27j. Lepton pairs are formed in a fire cylinder which undergoes 
a hydrodynamic expansion. The hadron-chemical analysis of a large body of hadronic 
heavy-ion data has shown that the final observed particle abundances at the SPS 
energies are consistent with a common chemical freezeout at temperatures of around 
175 MeV and baryon chemical potentials of around 270 MeV. In the two-freezeout sce- 
nario, during the subsequent expansion and cooling the system continues to strongly 
interact via elastic collisions maintaining the thermal equilibrium until the thermal 
freezeout. The absence of pion-number changing reactions then induces a finite pion 
chemical potential. To incorporate these features, Rapp et ai, have proposed a simple 
expanding fireball model I2Z] • This model assumes that the system is in thermal 
equilibrium up to time t max , when the thermal freeze-out occurs. The velocity field 
characterizing the expansion depends on space-time coordinates in the following way: 

v n (t,z) = (v z + a z t) — , (5.18) 

^max ("J 



v±(t,r) = (v r + a r t) — , 

^maxl^J 



where 



^max(i) = z + v z t + ^a z t 2 , (5.19) 

?"max(t) = r + V r t + ^a r t 2 , 

are the boundaries of the system at time t. We use the following parameters of 
expansion: t max = 11 fm is fireball lifetime, Zq = 4.55 fm is equivalent to formation 
time or initial temperature, r = 4.6 fm corresponds to the initial transverse nuclear 
overlap radius, v r = 0, and the parameters v z = 0.5, a z = 0.023 fm -1 and a r = 0.05 
fm -1 are adjusted to the final observed flow velocities. The time dependencies of the 
temperature and the baryon density are as follows 

T(t) = ( 21 M eV)ex P (- I gA_), (5 .20) 

260 . . 

= WY (5 - 21) 
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Figure 5.2: The assignment of momenta in the uo — > e + e~ decay, where g is the 
incoming momentum of the worp meson, qi and ^2 are the outgoing momenta of the 
leptons 



where 260 denotes the number of participating baryons, and 

V(t) = 2^ max (t)rL x (t), (5.22) 

is the volume of the fire-cylinder at time t. 

For the time dependence of the pion chemical potential, //(£), we assume, following 
Ref. jSll, the linear rise from 20 MeV at t = to 80 MeV at t = t max . 



5.4 Kinematics for dilepton yield 

Direct decays 

A dominant source of dileptons in the low-mass region is the direct leptonic decays 

of vector mesons. Although we focus on Dalitz decays for dileptons of mass in the 

region of interest, one also considers the p° and uj mesons, see Fig. 15.21 The dilepton 

rate formula for the p° — *> e + e~ decay is 

dN^- f d 3 q 1 3MT p 

-T p ^ e+e -f p (q ■ u{x))-— (5.23) 



(PxdAP J (2tt 3 ) "- e 6 JP ^ v 7r (rrip — M 2 ) 2 + (MT P ) 2 ' 
where M is the invariant mass of the lepton pair, x is a space time point, and T p is 
the total rho-meson width. The quantity r p _>/+/- is the decay width for p° — > e + e~ , 



- in ~ = 3 9p M» U " ' U ™ T (5 ' 24) 
where oiqed denotes the fine structure constant, and is the mass of the lepton. 
The function f p (q ■ u) in Eq. 15.231 is the thermal Bose-Einstein distribution of the p 
mesons in the thermal model, 

f p {q ■ u{x)) = exp[-( ^— )], (5.25) 
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Figure 5.3: The hydro dynamic expansion of the cylinder moving in the lab system 
with rapidity otpc m its own center-of-mass system. e + , e~ are the leptons formed in 
a fire cylinder, whereas vj_, vii are the velocities, which depend on space-time, and 
z max and r max are the boundaries of the system at time t; see section 4.3. 

with the pion chemical potential fi n . 
Expansion of the fireball 

To be as realistic as possible, we include the effects of the expansion of the medium 
formed in a relativistic heavy-ion collision, Fig 15.31 The lepton pairs are produced 
in a fire cylinder which moves in the lab system with rapidity olfc- The fire cylin- 
der expands hydrodynamically by itself. In particular, Eqs. 15.181 and 15.191 describe 
this expansion in the fire cylinder's center-of-mass system (CM). In such a case it 
is convenient to rewrite Eq. ()5.23|) in variables which are more convenient for the 
kinematics of the emission process and the geometry of the experiment. We introduce 
Mj_ = \J M 2 + q\ , the transverse mass of the dilepton pair, y lab , the rapidity of the 
pair measured in the lab system, uj_, the transverse four- velocity of the fluid element 
producing dileptons, and a lab , the rapidity of this fluid element in the lab. With these 
variables we have in the cylinder's rest frame 

q-u = M ± ^l + u 2 ± cosh (y lab - a lab ) - q ± • u ± . (5.26) 

The velocity of the fluid element in the lab is a relativistic superposition of the velocity 
of the fire cylinder in the lab and the hydrodynamic flow considered in the CM system. 



Thus we have 

iab i t>_i_cosh(a) 

a = a + ctFc = arctanh + ape-, u±= — . (5. 27) 

1 — vj_ cosh 2 (a) 

The velocities t>|| and t>j_ are now defined in the CM system and they depend on time 
and space coordinates. 

Next, we analyze the kinematic constraint of the CERES experiment, see Appendix 
C. The experimental acceptance cuts can be included using the function 

lab \ J d 2 qi±d 2 q2±d yi dy2 5{E q - E qi - E q2 )6®(q - qi - q 2 ) 



Acc(M,y iab ,q ± ) 



J d 2 q 1± d 2 q 2 ±dy 1 dy 2 5(E q - E qi - E q2 )5W(q - q x - q 2 ) 

(5.2? 



where qi2 are the momenta of the emitted electrons, yx,2 are the electron rapidities, 
and is a product of step functions which enforces the experimental setup conditions: 
2.1 =< yi j2 < 2.65, q^ 2 > 200 MeV, and 9 ee > 35mrad, where 9 ee is the angle between 
the directions of the leptons. We can assume that the rapidities and pseudorapidi- 
ties of the electrons are equal because of the smallness of the electron mass. The 
construction of the function Acc(M, y lab , q±) requires numerical calculation of a two- 
dimensional integral of a function involving a product of step functions, which can be 
done by a Monte Carlo method. With the inclusion of the experimental acceptance 
cuts, the dilepton production rate is 

dN e+e - f d 2 q ± 



d A xdM 



2Mjj^-J dy lab Acc(M,y la \q ± ) x 



1 3MT p 
X n(ml-M 2 y + (MT p y ^ l+l ~ tp (5 ' 29) 

To calculate the dilepton spectrum, we have to assume a model of the hydrody- 
namic expansion of the fire cylinder, as described in the previous subsection. Finally, 
the yield of dileptons produced during the expansion is 

T = / * / *"* / * (») • < 5 ' 30 ' 

-z max (t) 

where ^ d ^ dM ^ is given by Eq. ()5.29j) with elements defined in Eqs. I5.23j I5.25| 15.281 
Similarly, we calculate the dilepton rate formula for the ui — > l + l~ decay. 

Dalitz decays 

In order to calculate the Dalitz decay uj — ► 7r°e + e~, we choose the following kinematics, 
Fig 15.41 pi is the incoming momentum of the uj, q is the outgoing momentum of the 




Figure 5.4: The assignment of momenta in the uj — > 7r°e + e~ decay, where p\ is the 
incoming momentum of the u, p 2 is outgoing momentum of the pion, and q is the 
outgoing momentum of the virtual photon. The outgoing momenta of the leptons are 
indicated by q\ and q 2 . 



virtual photon, and P2 = P\ —q is the outgoing momentum of the pion, where q 2 = M 2 , 
and M denotes the invariant mass of the lepton pair. 

The particles are on the mass shell, hence p\ = m^, p\ = m 2 , and Pi ■ q = 
\{ m t + M 2 — ml). In Fig. 15 A[ q 1 and q 2 are the outgoing momenta of the leptons. 
The dilepton production rate from the Dalitz decays of the u is given by the general 
formula 

fjfj e + e - r r 

= E / dP MPi ■ u ) dp i t 1 + A(Pa • u )\ / dQ 1 dQ 2 (27r) 4 X 

ss' p 

5 {4) (pi-p 2 -gi-g 2 )|-M {p) | 2 (5.31) 

where x is a space-time point, fuj(pi ■ u) and /,r(p2 • u) are the distribution functions 
for the uj meson and the pion, respectively, u is the four velocity of the medium, 
dPi = (^Fpfg > d Qi = (2^)3^' ) l-^| 2 is the transition amplitude, ^ ss , is a sum over 
spins of the leptons, and J2 P denotes the sum over the polarizations of the decaying 
meson. The amplitude for these processes has the invariant form 

cm 1 

= A Jtr el ) (p 1 )G(q)—e-^eL ff (5.32) 
where A la is given by Eq. (j2.29|) . e 1 is the polarization vector of the u, G(q) is 



the denominator of the rho-meson propagator, is the conversion factor from 
the vector- dominance model, \ comes from the virtual photon propagator, and 
L a = u(qi, s) , ~f cr v(q2, s'), with u and v denoting the lepton spinors. We can write 
the amplitude in a more convenient notation 

MW = V®^L a (5.33) 
q z 



where 14 (p) is the ujttY vertex, denoted as V^> = V} p> G(q) e -^- , and V^ p> = A la eJ p) ( Pl ) 
is the ujp-K vertex. We obtain 



Hp) _ T/(p)^^^ e^ra £ 



\M<?>\ 2 = Vj p) V^ > -^L CT {L< 7 'y. 

Now, we rewrite the dilepton-rate formula Eq. 15.311 as 

dN e+e ~ 



d 4 



£ / dQ 1 dQ 2 (2n) 4 5^(q - 9l - g 2 )L CT (L CT ')*47ra e 



x 



(5.34) 



(5.35) 



£ / dP^Cp! • «)dP 2 [1 + /.(pa • u)] 5 (4) (Pi - pa - ff)V?Vj?>*i 
p ^ J 

where a e , m , = e 2 /An is the fine structure constant, L a (L a ')* depends on q± and g 2 , 
V^V^* depends on pi and p 2 , and M 2 = g 2 . First we compute the integral over 
dQi and in the above expression. It has the form 

S aa ' = -47TCW J dQ 1 dQ 2 (2ir) i S^(q - Ql - q 2 )L"(L*')* , (5.36) 

where explicitly, 

L\L°'y = u( gi , s) 7 CT v(g 2 , s')v(q 2 , s')Y'u( qi , s) = (5.37) 
Tr 



aQi+mi jqi-mi 
7 — ~ 7 



2m, 



2m, 



i tof + - 9 aa '(qi ■ <&)) - g c 



m 



In general, we can decompose the symmetric tensor S aa , which depends on the vector 
q, as 



= AQ aa + B 



q a q a 
^2- 



(5.38) 



where Q aa ' is a projection operator given by Eq. 12.571 However, due to gauge invari- 
ance, B = 0. Indeed, contracting S aa ' with q a q (7 'we find 



S^g^TTtw y rfg 1 dQ 2 (2 7 r) 4 5( 4 )(g - g x - g 2 ) 



mf 



\2q l - qq 2 - q- q 2 q l - q 2 ) - q 2 



Bq 2 



(5.39) 



The term in square brackets is due to the condition q = q± + q 2 and the lepton on 
the mass shell condition q\ = q\ = mf; hence, as required by general principles, 5 = 
. To extract A we contract S aa ' with g aa /, obtaining 

A = \s™'g aa , = 4ttcw / dQ 1 dQ 2 (27r) i 5^(q -q x - q 2 )-^(M 2 + 2m 2 ) (5.40) 
3 J 3mf 

The phase-space integral is most easily obtained in the rest frame of the decaying 
virtual photon of mass M. We get 



J dQ l dQ 2 (27r) i 5^(q - Ql - q 2 ) = ^(1 



■2 4 m 2 



M 2 



(5.41) 



and the partial decay width for the virtual photon is 
„ A a P „, ( 4m 2 



2M 3M V M 2 

Thus 

dN e+e ~ 



(^ + 2m 2 ). (5.42) 



d 4 x 



= £ / A / dPi/ w (pi • ^P 2 [1 + A(P2 • u)} 
p J J 

x 5^( Pl -p 2 -q)T p ( Pl ,q)^T r ^, (5.43) 



where we have defined 

T p ( Pl ,q) = V^V^*Q aa '- (5.44) 

A useful form of the dilepton-rate formula can be obtained through the use of the 
identity 

j d A q = j d A q j dM 2 5(M 2 - q 2 ) = J d 3 q J dM 2 ^. (5.45) 

We have 

dN e+e - f d 3 



d^xdM 2 

V 



■2 ■ U)\ X 



5 {4) (pi -P2- v)jp TP (p^ ?)r 7 ^ J+ i- = Yl [ d ® [ dp Mpi ■ u ) x 

v J J 

dP 2 [1 + U(p 2 ■ u)} (2n) 4 S^( Pl -p2- q)-±^T p ( Pu g)I>^ +J -. (5.46) 
Explicitly, the quantity T p ( P i, q) is equal to 

T p ( Pl: q) = A^^^O^v^S'bi)^^)^! 2 - (5-47) 

The transversely polarized uj has two helicity states, with projection equal to ±1 on 
the direction of q, while the longitudinally polarized u has one helicity state, with the 
corresponding projection equal to 0. Taking this into account we find: 



T ( P \, q) = L 11 Q aa Ay a A*, a ,, (5.48) 
and 

T T ( Pl , q) = T^Q™'A, a A;, al , (5.49) 

which we may rewrite as 

J2T p ( Pu q) = (T L + 2T T )\F\ 2 , (5.50) 



where 

9p 9p V( m p + M ) + M 2 TjJ 

and T p is the width of the p meson of virtual mass M. The width is dominated by 

2 3 

the p — > 7T7T channel, hence T p = (M 2 - 47712)2. Then we can write 



dN e+e f d 3 q f d 3 p\ £ { _ N f d 3 p 2 , > 4 



d 4 xdM 2 J (2n) 3 2E q J (2tt) 3 2E p 'J^ Pi ' ^ J (2vr) 3 2£ P2 ^ * 
[1 + U(p 2 ■ u)] 5{E P1 - E P2 - E q )6®(£ -n~~q)x 
{T L + 2T T )\F\ 2 -^T^ l+l -. (5.52) 

The distribution functions have the thermal form f u = f w {pi-u) = exp[— (jpx-u) /T(t)\ 
and / w = [1 + /^(pa .«)] = ! + exp[-(p 2 • u)/T(t)], where (p x • u) = E Pl u - Pipy - 



Pi±u±(cosa cos 9 — sin a sin 9) , and m = yl + + m| . Integrating our expression 
over the p 2 momentum and substituting d 3 g as q±dq±dq\\da and <i 3 pi aspij_dpij_dpi||d#, 
where the angle a is between q± and tTl, and between p 1 ± and g^, we obtain 

dN e + e - I" q ± dq±dq\\da f Pi±dpi±dpi\\d9 n 



d 4 xdM 2 J (2vr) 3 2£ g 7 (2vr) 3 2£ Pl JUJJn E p2 

1 

tt'M 3 



5(£ Pl - E P2 - E q )(T L + 2T T )\F\ 2 -±-T r ^ l+l -. (5.53) 



The on-mass-shell conditions are 



E P1 = Jml+ P 2 1± +p 2 w (5.54) 



E P2 = Jml +p 2 2 = Jml + (pi + q ) 2 = Wm 2 + (pi\\ - q\\) 2 + p 2 1± + q\ - 2p 1± q ± cos9 



E q = ^M 2 + g|+g 2 . 

Thanks to energy conservation we can eliminate integration over the angle 9, from 
the Dirac delta in Eq. 15.531 We get 5(9 — 9 ), where 



-(M 2 + m 2 - ml + 2p 1{m - 2Jml+p 2 + p 2 /M 2 + q 2 + q\) 

cos9 Q = 1 ¥ . (5.55) 

2pi±g± 

We have to remember that cos# must be in the range [—1,1]. For dilepton 
production from tu-meson decays we can write: 

dN e+e ~ f Trdq ± dq ] \da dp 1± dp Ml 2 6(1- | cosflp |) 1 2 
d 4 xdM 2 J {2ix)nE q {2ix)nE pi Tuh |sin^ | irM 3 L 1 

(5.56) 

where the factor 2 in the numerator comes from the two solutions of the equation 
cos# = COs6> . 
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Figure 5.5: The dilepton yield for the direct decays u> — > e + e~ and p — > e + e~, and the 
Dalitz decays u — > 7r°e + e~ and p a — > n a e + e~ , plotted as a function of the invariant 
mass, M e+e _. The two upper plots are for the vacuum value of the coupling constant; 
lower plots include the medium-modified irup vertex. We show our results for two 
different temperatures, T = 210 MeV (left), from Ref. [H] and T = 150 MeV (right), 
|l()9j . The thick line indicates the sum of all contributions. 



Next, we apply the experimental acceptance cuts, which are included in the func- 



tion 



Acc{M, q h q ± ) 



j d 2 q l± d 2 q 2 ±dqi\\dq 2 \\^E-; S(E q - E qi - E q2 )8^(q - qi - q 2 ) 
/ d 2 q 1± d 2 q 2 ±dqi\\dq 2 \\ 



E ql E q2 



5(E q -E qi -E q2 )5^(q-q 1 -q 2 ) 



(5.57) 



where the electron rapidities, see eq. 15.281 were rewritten as dy = dqn/E q , with 
E q = \J M 2 + q\ + q 2 . Finally, after using the fire-cylinder hydrodynamic expansion 
model of [211 |22j described in sec. 4.3, the yield of dileptons produced during the 
expansion is (the Bose enhancement is neglected for brevity) 

dN S ~^ e ~ /''max /TmaxW 

" = Jo 



dt 



27frdr 1 dz \ M ^n,7 AcC ( M > % > 9-J-. 

-max 

(*) 



dM 2 J J \d A xdM 2 ' 
In the case of the Dalitz decay of the p meson we apply the same formulas. 



(5.5* 



Chapter 6 

Results of our model 



In the following we present our numerical results. We begin with calculations without 
expansion of the fire-cylinder and without the acceptance cuts, Eqs. 15.281 and I5.57| 
shown in Fig. 15.51 We do this in case to compare with other works. New and impor- 
tant result of this study is to applying in our calculations the CERES experimental 
kinematic cuts, what in other theoretical works is usually skipped. 

In Fig. 15.51 we plot the dilepton yield as a function of the invariant mass, M e+e _. 
On the left-hand side we show the results for the temperature of 210 MeV of Ref. 
|2E]; see section 4.3. The upper plot is for the vacuum value of the nup coupling 
constant. We display the direct contributions, such as uj — > e + e~ and p — > e + e~, 
which are indicated by pink and blue curves, and the Dalitz contributions of both uj 
(black curve) and p (green curve) mesons. On the lower plot we present the same 
results but with medium modified irujp vertex. The thick red line indicates the sum 
of all direct and Dalitz contributions. We see that the Dalitz decays are enhanced 
by about a factor of 2 compared to the vacuum values. On the right-hand side we 
show the results for a lower temperature of 150 MeV, with similar conclusions. It is 
visible that the direct decays are dominant in the region above 0.35 GeV and that 
the number of produced particles increases with growing temperature. Here we also 
observe the peak from the direct decay of the omega meson, which is not observed in 
the data because of the too low experimental resolution. 

Next, we include in our calculations the expansion effects and the CERES ac- 
ceptance function. In Fig. 16.11 we illustrate the dilepton yield for the uj — > 7r°e + e~ 
(left-hand side) and p a — > 7i a e + e~ (right-hand side) decays, plotted as a function of 
the invariant mass of dileptons, M e+e _. The plots are done for the maximum fireball 
temperature of 210 MeV and for different values of the baryon density, ps(t) (green 
dotted curve) and 2 ps(i) (red curve). In our calculations we use the Rapp-Wambach 
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Figure 6.1: The dilepton yield for the Dalitz decays. The dashed lines are for the vac- 
uum value of coupling constant, whereas the solid lines include the medium modified 
irujp vertex. Green and red lines show the results for different value of baryon density 
Psit) and 2 psit), respectively. We show our results for T = 210 MeV. 

model of the fireball, but we change the temperature from T=210 to T=300 and in- 
crease two times the baryon density. This significantly enhances the medium effects. 
By the blue dashed line we indicate the case with the vacuum value of the Troop vertex. 
We can observe that the medium effects from the Dalitz decays of vector mesons are 
clearly visible, over the entire region of M e+e _ = 0.2 — 0.6 GeV. For both Dalitz 
decays u — > 7t°e + e~~ and p a — > n a e + e~ there is a difference of about a factor of 2 for 
the case with normal baryon density compared to the vacuum value. Simultaneously, 
we notice that this enhancement factor rises to a higher value of about a factor of 5 
with the growing fireball density. 

We have a similar situation in Fig. 16.21 where the same results as in Fig. 16.11 but 
for temperatures of 300 MeV, are presented. The dilepton production from Dalitz 
decays is found to be strongly enhanced for the medium-modified vertex compared 
to the vacuum for this value of the temperature. This difference grows even larger, 
to about one order of magnitude, with the growing baryon density. It is naturally 
explained since we expect to observe more particles produced at higher temperatures 
than at lower ones. 

We continue our discussion with Fig. 16.31 and 16.41 where we present the dilepton 
spectra from the direct and Dalitz decays compared to the CERES experimental data, 
taken from Ref. [EH], for two different temperatures, T=210 MeV and T=300 MeV. 
Both plots are for the same density equal to twice the normal baryon density. As we 
can see the medium modifications help, but not enough to reproduce the experimental 
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Figure 6.2: Same as Fig. 16.11 but for the T = 300 MeV. Left side, the dilepton yield 
for the Dalitz decay uj — > 7i°e + e~; right side, for p a — > n a e + e~ , plotted as a function 
of the, M e+e _. Dashed lines correspond to the vacuum value of the coupling constant, 
solid lines include medium modified nujp vertex. Green and red lines show the results 
for different values of the baryon density psif) and 2 respectively. 

data. However, for the higher temperature of the fireball, T=300 MeV, results lie 
much closer to the data. We notice that the medium effects from Dalitz decays of the 
uj and p mesons are large and are specifically about two times larger than in vacuum. 
We see that the contributions from the p a — > ir a e + e~ decay are equally important as 
the ones from the Dalitz decay of the uj meson and should not be neglected during the 
considerations connected with the low-mass dilepton spectrum, as is frequently done. 
It is visible from our analysis that in the region of interest, 0.2 — 0.6 GeV, where the 
existing calculations have problems explaining the experimental data, our effect helps 
but still fails to reproduce the experimental data. 

We have to admit that there are other processes which contribute in the explored 
region of M e+e _, like Dalitz decays of rj and rj mesons which we do not include in 
our calculation. We have not taken into account other effects such as shift of the p 
meson peak position or broadening of the widths. We also do not show the 'cocktail' 
contribution of decays of hadrons on top of our results, which would not be consistent. 

After all, we stress that the nujp coupling constant undergo substantial modifi- 
cation in the nuclear medium and directly influence the dilepton production rate, to 
which this thesis was devoted. 
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Figure 6.3: Dilepton emission rate for 158 AGeV Pb+Au CERES experiment from 
the Dalitz and direct decays (for p± > 0.2 GeV). The dashed and solid line correspond 
to the u — > 7r°e + e~ and p a — > 7r a e + e~ decays, respectively. The solid blue and green 
line correspond to direct contributions p — > e + e~ and — > e + e~ . We show our results 
for temperature T = 210 MeV and double normal value of the baryon density. 
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Figure 6.4: Same as Fig. O but for T = 300 MeV. 



Chapter 7 

Summary and discussion 



In this thesis we have considered the modifications of the ncop coupling in dense 
nuclear matter and its influence on dilepton spectra in relativistic heavy-ion collisions 
coming from the Dalitz decays. 

In the first part we have analyzed the noup vertex in the context of the Dalitz de- 
cays of 7T°, uj and p particles. Our calculations were done in the framework of a fully 
relativistic hadronic theory with mesons, nucleons and A isobars. We have presented 
diagrams, propagators, typical vertices for the interactions with and without A res- 
onance, and the in-medium tensor structure of the nup coupling constant. Presence 
of the additional four-vector, the four-velocity of the medium, made this structure 
more complicated. We have worked at zero temperature and in the leading-density 
approximation. We have considered decays of particles in the rest and the moving 
frame with respect to the medium. 

In the case where the decaying particles are moving with respect to the medium 
we have analyzed the kinematics in a convenient way, in the rest frame of the nuclear 
matter, not in the rest frame of decaying particle, as it frequently done in the vacuum 
calculations. Because each polarization of the decaying particles behaves differently 
the properties of these particles are different, in particular their widths. Thus, we have 
analyzed the dependence of the widths for transverse and longitudinal polarizations 
on the invariant mass of dileptons. 

The main conclusion of our investigation is that the medium effects on the nup 
coupling are large and primarily come from the process where the A is excited in 
the intermediate state. On the other hand, for the ir° — > 77* decay the effects of 
the A isobar cancel almost exactly the nucleon particle-hole excitation, such that 
the medium contribution is small. The dilepton yield from the ir° decays are not 
changed by the medium, because virtually all pions, due to their long lifetime, decay 
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outside the fireball and because the width for this decay is practically unaltered by 
the nuclear matter. However, we have found a sizeable increase of the nup vertex 
compared to its vacuum value for the Dalitz decays from u and p mesons. This 
coupling constant is enhanced by about a factor of 2 for the u — > 7r°7* decay and by 
a factor of 5 for the p a — > 7r a 7* decay. For the decaying particle moving with respect 
to the nuclear matter, the medium effects decrease with growing magnitude of three 
momentum, but remain significant for a momenta up to ~ 200 MeV, the relevant 
value for temperatures typical in heavy-ion collisions. This results confirm that the 
medium effects are substantial. 

Through-out our calculations we were following a relativistic covariant approach. 
We have also shown in detail how to obtain the leading-density approximation for 
dileptons involving loops with density-dependent nucleon propagators. In order to 
simplify our approach we have not included any form factors in the vertices, because 
they lead to fundamental problems with gauge symmetry conservation, in particu- 
lar the Ward-Takahashi identities and current conservation are violated. We have 
considered the Dalitz decay of the p meson which was never analyzed before. 

In the second part we have applied our model to evaluate the dilepton production 
rate from the Dalitz decays of vector mesons in relativistic heavy ion collisions. In 
order to estimate the dilepton yields we have used the VDM and including expan- 
sion of the medium. In our calculations we have adopted the Rapp-Wambach model 
of the hydrodynamic expansion of the fire-cylinder, which includes longitudinal and 
transverse expansion. The kinematic constraints of the CERES experiment also is 
used. The construction of the acceptance function required a numerical calculation 
of a two-dimensional integral involving the product of step functions, which we have 
obtained by a Monte Carlo method. 

In order to compare our results with other works, we have shown the dilepton 
yield as a function of the invariant mass for calculations without expansion of the 
fire-cylinder and without the acceptance cuts. We have presented plots for two temper- 
atures, 210 MeV and 150 MeV, with conclusions that the Dalitz decays are enhanced 
by about a factor of 2 compared to the vacuum and that the number of particles 
increases with temperature. On the other hand, the direct decays are dominant in 
the region around 0.4 GeV. 

Next, we have included the expansion and the CERES acceptance function. Our 
final plots are done for two fireball temperatures (T = 210 MeV and T = 300 MeV) 
and for density equal to twice of the normal baryon density. Comparing our results 
to the dilepton invariant mass spectrum measured by the CERES collaboration at 



CERN SPS we have found that the Dalitz decays are enhanced in the region of 
0.2 — 0.6 GeV, and medium modifications of coupling constant increase the effect. 
It is worth pointing out that the contribution from Dalitz decay of the p meson is 
equally important as the u>. However our results still underestimate the experimental 
data. 

Our model together with other effects, which are not considered here, such as 
creation of dileptons from other decays (77 and 77'), dropping masses or broadening 
widths, may be able to explain the present experimental data. The other effects 
should be taken into account in future work. 

For our numerical analysis we have developed a program written in 'Mathematica' 
with the very useful FeynCalc package for algebraic calculation in elementary particle 
physics. We have included in the Mathematica code the acceptance function, which 
was calculated in Fortran. 

The low-mass dilepton research is a quickly evaluating branch of theoretical and 
experimental physics. Inclusion of more sources of dileptons or detailed verifications 
of existed theoretical models, together with new experimental data will provide new 
insight into that relevant problem in relativistic heavy-ion collisions physics. 



Chapter 8 
Appendices 



8.1 Appendix A 

The contributions to B from diagrams (a), (b) and (c) of Fig. 12. HI for the process 
7T° — ► 7(^)7* are as follows: 



l7 r» m N (K p + l)gA9 P 9u,ml(3x 2 + m|) 

(4m^ — m^.)((m7vx 2 + m^m 2 ) 2 — x A m 2 ) ' 

W ~ gMimTvflm^x 2 + m N ml) 2 - (m 2 N - M 2 + x 2 ) 2 ™ 2 } [[mN a) * 

(m 2 ^ + m N M A - M A + x 2 )ml - (m N x 2 (m 3 N - m N M A M A x 2 )) - m 2 N ml), 



.*» 9nA7t9na p 9^F w (-(m 2 (Om + 

W 27Mim n [{m 2 N -Ml + x 2 ) 2 ml-{m N ml + m N x 2 ) 2 } y 1 M ^ 

+3M A )x 6 ) + +x 2 (4m^ + 6m^M A - 8m 3 N x 2 + 3m^(M^ - 3M A x 2 ) + 

+2m N (AM A - 2M\x 2 + x 4 ) + 3M A (M A - 3M A x 2 + 2x 4 ))m 2 + (4m^ + 

+6m%M A + 11M A - 17M A x 2 + 6M A x 4 - 2m^(4M A + x 2 ) - tt&(5M£ + 

+9M A x 2 ) + 4mjv(4M A - 5M A x 2 + 2x 4 ))m 4 n - m N (Am 2 N + ?>m^M A + 

+2x 2 )m 6 J), 

(8.1) 



where x = 
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The contributions to B from 7r° — ► 7(p)7* decay are as follows: 

B 7T°( P ) gAgpgojfnl(2m 2 N x 2 (3 + k p ) + [x 2 k p + 2m 2 N (l + K p ))ml) 

w 2mjv(4m^ — m 2 )((mj^x 2 + m^m 2 ) 2 — x^m 2 ) 

B ^{ P ) gNAngNApgujF^m^jx 2 - ml)(m 2 N x A (m 3 N - 2Mx 2 - m N {M 2 A + 2x 2 )) + x 2 

"'' 9^a((^a — m N) 2m l ~~ (itinx 2 — mNm 2 ) 2 )({rriNX 2 + m^m 2 ^) 2 — x 4 m 2 ) 

{Am% + m%M A + 3m N M A + 2M 5 A + m%(-7M A + 6x 2 ) - 3m 2 N (M 3 A - 

-M A x 2 ))ml ~ m N(m% + m 3 N M A + M A + m 2 N (-2M 2 A + 5x 2 ) + 

+m N {-M A + M A x 2 ))mi + m%m 6 n ), 
B n°( P ) gNAngNA P gujF v (-(m 2 N (2m N + 3M A )x 6 ) + x 2 (Am% + 6m%M + 8m N M A 
( c > 27M A m v ((M A — m 2 N ) 2 m 2 — {m N x 2 — m N m 2 ) 2 ) 

+3M 5 A + 8m 3 N x 2 + +3m 2 N (M 3 A + 3M A x 2 ))ml - (Am% + 6m%M A + 

+ 16m N M A + 11M A - 2m 3 N (AM 2 A - 5x 2 ) + m 2 N (-hM 3 A + 9M A x 2 ))m 4 n + 

+m 2 N (4m N + 3M A )ml) 

(8.2) 

with x = m p . 



8.2 Appendix B 

Numerators and denominators of Eq. 12.331 

For the case q = the contributions to B from Eq|23Sl corresponding to the diagrams 
(a,b,c) from Fig l2.3| can be written as 

(8 ' 3) 

The formulas for Nu\ and Du\ are very long. Numerators Nr^ become manageable 
in the formal case T A = 0. Below we show the numerators for the diagrams with 
nucleons only, because those with A isobar are too long. We find 



Nfa = gAg P guj(m 2 N K p m n 8 (m 2 N - mj) + 2m 2 N mJ s {-2m 2 N x 2 k p + (x 2 K p - 3m 2 N x 
x(l + Kp^mJ 2 + KpTriu*) - m 2 N (x 2 - mj 2 f \-(m 2 N x A k p ) + x 2 {x 2 k p + 
+m 2 N (Q + AK p ))m w 2 + m 2 N (2 + K p )m UJ A ) + m^(6m%x 4 K p + 2m 2 N x 2 (-(x 2 K p ) + 
+3m 2 N (l + Kp^mJ 2 + (2m 2 N x 2 + x a k p + 2m%(5 + ^K p ))mJ t - m 2 N K p mJ > ) - 
—2m 2 N rn 1T 2 (2rn 2 N x e 'k p — x (x 2 n p + 3m 2 N {l + K p ))m u 2 + x 2 {— x 2 + 2m 2 N (5 + 
+AK p ))mJ + (m 2 N - x 2 )(l + K p )mJ)), 



N ( a ) = -(gAg P g^m p 2 (6m%(l + n p )mj + m%(x 2 -m p 2 ) (2m 2 N x 2 (3 + k p ) + (x 2 k p + 
+2m 2 N (l + K p ))m 2 ) - m^{2m%x 2 {3 + 5k p ) + m 2 N (x 2 (2 - 3k p ) + 10m 2 N x 
x (1 + K p ))m 2 + x 2 K p m p 4 ) + 2m 2 N m 7T 2 (m 2 N x A (-3 + k p ) - x 2 (x 2 (l + 
+2k p ) - 2m 2 N (5 + 3K p ))m 2 + (-x 2 + m 2 N (l + K p ))m A ))). (8.4) 
The expressions for Di for the case r A = for u — > 7r°7* decay, where x = m p , are 
Dfy = 2(2m N - m u) )(2m N + m u) )(m N x 2 - m N m 2 - x 2 m w + m N mJ 2 )(m N x 2 - 
-mNin-K 2 + x 2 m ul + mNm IJ} 2 )(m^ 2 {mN - m u ) + m N (—x 2 + m w 2 )) x 
x(m n 2 (m N + mj + m N (-x 2 + mj 2 )), 

Ufa = 9M 2 m^(-(m N x 2 ) + m N m w 2 + (m 2 N - M A + x 2 )m^ - m N m^ 2 )(m N x 2 - 

—rriNm^ + [m N — M + x )m u) + m^m^ ){m^x + (m^, — M A )m w — mjv x 
xm^ 2 + m n 2 (-m N + m aJ ))(-(m Ar a; 2 ) + (m^ - M^)m w + m N m^ 2 + 
+m 7T 2 (m N + mo,), 

D^,) = 27 M^m^(-(m N x 2 ) + m N m w 2 + (m 2 N - + x 2 )m^ - m N m^ 2 )(m N x 2 - 
-mum^ 2 + (m^ - M A + £ 2 )m w + m N m UJ 2 )(m N x 2 + (m 2 N - M 2 A )m w - 
-rriNm^ 2 + m n 2 (—m N + m w ))(-(TO;v£ 2 ) + (m 2 ^ - M^)m w + m^m^ 2 + 
+m w 2 (m N + mj) 

(8.5) 

The expressions for Dj for the case r A = are for the p — > 7r°7* decay, where x = 



= 2(2m N — m p ){2m N + m p ){m N x 2 — m N m 7T 2 — x 2 m p + m N m 2 ){m N x 2 — 
—m^m^ 2 + x 2 m p + m^m^im^im^ — m p ) + m^i—x 2 + m 2 )) x 
x (m 7r 2 (m A r + ra p ) + m N (-x 2 + m p 2 )), 

Z?^ = 9M A (m N + M A — m p )(m N + M A + m p )(m N x 2 — m N m w 2 — x 2 m p + 

+mNiri p 2 )(m]yx 2 — m^m^ 2 + x 2 m p + m^m 2 ){m^x 2 + (m 2 ^ — M A )m p — 
-m N m 2 + m 2 {—m N + m p ))(—(m N x 2 ) + (m^ - M A )m p + m N m 2 + 
+m n 2 (m N + m p )), 

Z^ c) = 27M A (m N + M A - m p )m p (m N + M A + m p )(m N x 2 + (m 2 N - M A )m p - 
-m N m 2 + m 7r 2 (-m A r + m p ))(-(m N x 2 ) + (m 2 ^ - M A )m p + m N m 2 + 
+m 7r 2 (m A r + m p )). 

(8.6) 



8.3 Appendix C 



Ceres Acceptance function 

In this appendix we develop the expression for the acceptance function of the CERES 
experiment. The kinematics cuts include the pseudo rapidity range, reflecting the 
angular acceptance of the detector cuts for transverse momenta of the measured e + 
and e~, as well as the minimum angle between the directions of the e + and e~ . The 
kinematics of the decay is depicted in Fig. 18.11 




Figure 8.1: The decay of 7* into a lepton pair. 



We use the energy and momentum conservation 



where 



E = Ei + E2, 

Q\\ = P||i+P||2, (8.7) 



E = m_LCOsh(y), Ei >2 = p±i,2C0sh(y 1>2 ) 
q\\ = m_Lsinh(y), p|| 1(2 = p±i, 2 sinh(y 1)2 ) 



with y denoting the rapidity of 7*, and m± = a/ m 2 + q\ . Next, from the relation q± 
P±i + P±2 w e obtain 



RL2 = y q\ + P±i 2 - 2q±P±i cos ^, (8.9) 

where £ is the angle between q± and p±i. Solving Eqs. 18.71 with the help of Eq. 
we obtain the solution for the rapidities of e + and e~ of the form 

sinh(?/ 1 ) (1) ' (2) = - — smh(y)(m 2 ± + p 2 ±1 - p\ 2 ) ± 

2m_Lp_Li 



± Jmi-2(pl 1 +pl 2 )ml + (pl 1 -pl 2 ycosh(y), (8.10) 



2.6° 



Figure 8.2: The Ceres acceptance function for different value of M, left-hand side 
M = 0.2, right-hand side M = 0.6. 



sinh(y 2 ) (1) ' (2) = ~ — smh(y)(m 2 ± - p 2 ±1 + p 2 ±2 ) ± 

-in .P.i 



and 



± yjm 4 ± - 2{p 2 ±1 + p 2 ±2 )m 2 ± + (p 2 ±1 -p 2 ±2 ) 2 cosh(y). (8.11) 

The (1) and (2) label the two possible branches (the situation is similar to the case 
discussed in Chapter 2, section 2.3.) These solutions are real if the condition 



'mi - 2{p\ l + pl 2 )m\ + (p 2 ±1 - p\ 2 f > 0, (8.12) 

is fulfilled. The upper limit for the integration over p±i we obtain when £ = in 
Eq.EH Then 

P±i = \(q± + \Jm + ql). (8.13) 

The Ceres acceptance function is 

/ d 2 p l± d 2 p 2 ±dy 1 dy 2 4> 5{E q - E Pl - E P2 )5 {3) {q - Pi - p 2 ) 



Acc(M, y, q ± ) 



2tt 

(8.14) 



Using 5 (3) (q — pi — p 2 ) = <5 (2) (q± - p±i - p_L2)<$(qj| _ P||i _ Pp), where the two 
dimensional delta eliminates integration over p\ 2l we are left with 

d(E q - E pi - E P2 )d(qii - pp - pp) - — ; r . (8.15) 

p±iP±2COsh(y 1 - y 2 ) 

In Eq. I8.14[ is the product of step functions which includes the experimental cuts, 
such as 2.1 < yi t 2 < 2.65, p±i j2 > 200 MeV, 8 CC > 35 mrad. It has the explicit form 



x 6 



6(3/1 - 2.1)0(2.65 - yx)e(y 2 - 2.1)(2.65 - y 2 )&( P±1 - 0.2)9(p ±2 - 0.2) x 
RliP_L2 +p±isinhyip_ L2 sinh?/ 2 



cos (0.35) 



RLicoshyipj_2Coshy 2 



8.16) 



As an example in Fig. 18. 2[ we show a plot of Acc(M, y, q±) for different values of mass 
M. We calculate the CERES acceptance function using the Monte Carlo integration 
method in Fortran. 



8.4 Appendix D 

We are going to proof that the diagrams of Fig. l2.3l with more than one Sd propagator 
vanish for kinematic reasons. We write the two nucleon density propagators, Sm and 
Sd2 in the following way 

Sdi = 1(7% + m N )—8(k - E k )9{k F - |k|)], (8.17) 



27T 

S D 2 = 2(7% + 7% + m N )- — 5(ko + p - E k+p )6(k F - |k + p|)]. (8.18) 
From above expressions 18.171 and 18.181 we have 



k = E k = ^/k 2 + m% and k + p = E k+P = yj (k + p) 2 + m%, (8.19) 
thus, after simply calculations, we obtain 

m\ + 2p h = 2 k p cos 6, (8.20) 
with p = E p = ^/p 2 + m 2 . We have that 



ml + 2 A /p 2 + m 2 \/k 2 + mL 
cos 9 = -Z ov > 1, (8-21) 



which ends the proof. 
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